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ABSTRACT 

Different theoretical methodologies lead to order-of-magnitude variations in predicted galaxy-galaxy merger 
rates. We examine how this arises, and quantify the dominant uncertainties. Modeling of dark matter and galaxy 
inspiral/merger times contribute factor ^ 2 uncertainties. Different estimates of the halo-halo merger rate, the 
subhalo "destruction" rate, and the halo merger rate with some dynamical friction time delay for galaxy-galaxy 
mergers, agree to within this factor ^ 2, provided proper care is taken to define mergers consistently. There are 
some caveats: if halo/subhalo masses are not appropriately defined the major merger rate can be dramatically 
suppressed, and in models with "orphan" galaxies and under-resolved subhalos the merger timescale can be 
severely over-estimated. The dominant differences in galaxy-galaxy merger rates between models owe to the 
treatment of the baryonic physics. Cosmological hydrodynamic simulations without strong feedback and some 
older semi-analytic models, with known discrepancies in mass functions, can be biased by large factors 5) 
in predicted merger rates. However, provided that models yield a reasonable match to the total galaxy mass 
function, the differences in properties of central galaxies are sufficiently small to alone contribute small (factor 
^ 1.5) additional systematics to merger rate predictions. But variations in the baryonic physics of satellite 
galaxies in models can also have a dramatic effect on merger rates. The well-known problem of satellite 
"over-quenching" in most current semi-analytic models (SAMs) - whereby SAM satellite populations are too 
efficiently stripped of their gas - could lead to order of magnitude under-estimates of merger rates for low- 
mass, gas-rich galaxies. Models in which the masses of satellites are fixed by observations (or SAMs adjusted 
to resolve this "over-quenching") tend to predict higher merger rates, but with factor ^ 2 uncertainties stemming 
from the uncertainty in those observations. The choice of mass used to defined "major" and "minor" mergers 
also matters: stellar-stellar major mergers can be more or less abundant than halo-halo major mergers by an 
order of magnitude. At low masses, most true major mergers (mass ratio defined in terms of their baryonic or 
dynamical mass) will appear to be minor mergers in their stellar mass ratio - observations and models using 
just stellar criteria could underestimate major merger rates by factors ~ 3 — 5. We discuss the uncertainties in 
relating any merger rate to spheroid formation (in observations or theory): in order to achieve better than factor 
^ 3 accuracy, it is necessary to account for the distribution of merger orbital parameters, gas fractions, and the 
full efficiency of merger-induced effects as a function of mass ratio. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: active — cosmology: theory 



1. INTRODUCTION 

In the now established ACDM cosmology, structure grows 
hierarchically (e.g. White & Rees 1978), making mergers an 
inescapable element in galaxy formation. The galaxy-galaxy 
merger rate, as a function of properties such as galaxy mass, 
redshift, and mass ratio, is a quantity of fundamental inter- 
est. It is critical for informing models of the growth and 
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assembly of galaxies, the distribution of bulge and spheroid 
mass, properties of disks such as their thickness, morpholo- 
gies, and flaring, and the fueling and growth of the most lumi- 
nous starbursts and infrared systems as well as massive BHs 
and quasars. 

The last few years have seen the emergenc e of a "concor- 
dance " precision ACDM cosmology (see e.g. iKomatsu et al.l 
120091 and refer ences therein), with remarkable convergence 
between different probes of structure formation. Mean- 
while, cosmological dark-matter simulations have developed 
the ability to track large populations of dark matter halos 
and substructure within halos over cosmological volumes and 
timescales. Different calculations now yield results for e.g. 
the dark matt er halo mass funct ion that agree at the ~ 5 — 10% 
level (see e.g. lReed et al.ll2007i and references therein). 

Despite these advances, however, there has yet to be a 
similar convergence in theoretical predictions of the galaxy 
merger rate. Depending on the models used to convert dark 
matter halo growth histories into galaxy growth histories, dif- 
ferent theoretical predictors of the merger rate in the same 
cosmology can differ by more than an order of magnitude (see 
e.g. the comparison in Jogee et al. 2008; Bertone & Conselic^ 
120091: iLopez-Sanjuan et aLll2009al) . And the differences can- 
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not be trivially attributed to e.g. the known fact that models 
predict somewhat different absolute galaxy masses and abun- 
dances - in terms of mergers per galaxy, there is similar vari- 
ation. 

Some of this may owe to the fact that dark matter merger 
rates are - unlike dark matter mass functions - derivative (in- 
stantaneous) quantities and so are more sensitive to choices of 
definition and methodology in simulations. But much of the 
difference owes to the fact that the mapping from dark matter 
merger rates to galaxy merger rates is non-trivial, and requires 
some relation between the history of a given dark matter halo 
and a galaxy inside that halo. In models that attempt to pre- 
dict galaxy formation in an a priori sense, from physical pre- 
scriptions for gas cooling, star formation, and feedback from 
stars and accreting BHs, the resulting predicted halo to galaxy 
"map" and corresponding merger rates can be very sensitive 
to the input prescriptions. It has been shown that even in 
otherwise identical models, small changes in the treatment 
of gas cooling or stellar wind physics can lead to order-of- 
magnitude changes in the predicted major merger rates : com- 
pare th e merger rates in Bower e t al. (2006) and Font et all 
hoOSi) (our §l433Tl. or those in IdeXu cia & Blaizot (2001 
and I fiertone et alJ JioOTh (presented in iBertone & Consehca 
120091k 

These quantitative differences in the number of mergers per 
galaxy have led to qualitatively different conclusions regard- 
ing galaxy formation. Most models, especially those based on 
empirical halo occupation constraints or direct cosmological 
hydrodynamic simulation, have found that there are too many 
mergers in low mass systems to explain the survival and 
preva lence of disks in the simplest scenario (Granatoet al] 
20041: ISom erville et al. 2001, 2008a; Kodaetal. 2009; 



Stewart et a l. 2009a; Khochfar & Silk 2009 ; .Hopkins et ah 
2009bl; [Ste wart et al. 2009b; Sommer-Lar sen et alj 120031; 
Abadi et al. 20 03t iGovernato et a l. 2007; Robertson et al.' 
2Q04; Okam otoetalJ l2005t "jSc annapieco et al. 2008; 
Kazantzidis et ani2008ll2009tlPurcell et alJl20091) . 



More directly, observed merger rates find that ^ 5 — 
10% of low or intermediate-mass galaxies (M* < IO'^Mq 
at redshifts z ^ 0.2 — 1.2) are in mergers - strongly mor- 
phologically disturbed or in "major" similar-mass pairs 
about to merge (at small scales and small relative ve- 



locities) ([Bridge et al.l l2007l 120 lO t 'Karta ltepe etH 



2008b 



ConseUce etal.) 120091: lJogeeet al.1 2009; L otz et al.1 
Lin et all 120081 iRobaina et al.l 1201 0). This is equal to or 
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higher than the predicted rates from the models above. 
Even without correcting for the short expecte d duty cy- 
cles/lifetimes of this merge r phase (see e.g. iLotz et aH 
l2008al: IConsehce et al.l l2009l) . these fractions are already 
higher than the observed fraction of bulge-dominated galax- 
ies at these m asses (Allen et al. 2006; Benson et al. 2007; 
iDominguez ^lmero et al.ll200 8; Weinz iri et al.l l2009). Like- 
wise the number density of IR luminous systems appar- 
ently in mergers at high redshifts is so high as to yield 
tension with the relic bulge mass in the Universe (see 
e.g. ^ unger et al. 2007; Genel et al. 20081; jTacconi et alj 
l2008rTHopkins et al.ll2006t ICimatti et al . 2008). Typically, 
some additional physics (often some form of strong stel- 
lar feedback) must be invoked to preserve high gas frac- 
tions and suppress the efficiency of bulge formation. In- 
deed, this is always the case in hydrodynamic simulations 
of galaxy formation. This has led to the developrnent of 
various feedback m odels (e.g. | Sommer- Larsen et al.1 120031: 
lAbadi et al.1 l2003t IGovernato et all 120071; lOkamoto et alJ 



120051: IScannapieco et al.l l2008l and references therein) 
and emphasis on gas-richness as a stab i lizing factor 
("Robertson et al.l '2006^, 'Hopk ins et al] l2009al; IStewart et al] 
2009b; Moster et al. 2010a) that allows disks to survive such 
mergers. However, some semi-analytic calculations have 
claimed the opposite: that there are not enough mergers at 
low mas ses to explain eve n the small bulges observed at these 
masses dParrv et al.ll2009l) . 

On the other hand, at > L*, most empirical and semi- 
analytic models predict similar numbers of mergers and agree 
that ma ssive bulges are primar ily formed in major mergers 
(see e .g. lHopkins et al . 2010; P arry et al.l2009HCattaneo et all 
1201 Ol and references therein). But some simulations have 
claimed that there ar e more minor an d fewer major mergers 
at these high masses iNaab et al.ll2007l) . 

As observations of e.g. pair and morphologically dis- 
turbed fractions improve, it is rapidly becoming possible 
to empirically constrain merger rates and fractions at a 
level better than the apparent order-of-mag nitude scatter in 
predictions (see e.g. Kartaltepe et al. 2007; Lin et al. 20081 
Dpmingue et al.l |2009; Bundv et al. 2009; Conselice et ap 
2009; Robainaetal. 2009, 2010; Comerford et al. 200 J 
Lotz et al.. .2008h: .Bridge et al.i i2007j ; Lopez-Sanjuan et al] 
2009bl) . However, the conversion of such observations 



to merger rates still requires knowledge of appropriate 
timescales and contamination from projection effects and/or 
minor mergers - all of which mean that interpretation of the 
observations is somewhat dependent on how well or poorly 
theoretical predictions agree on basic quantities. Moreover, 
an observed pair or morphologically disturbed fraction is fun- 
damentally not invertible at much better than the factor of ^ 2 
level, as disturbances and dynamics are non-unique and de- 
pend on other quantities such as galaxy gas content; as such, 
they must be forward-modeled. 

Fundamentally, the galaxy-galaxy merger rate depends on 
two basic quantities: the dark matter halo merger rate and 
the manner in which galaxies populate halos. However, mod- 
els arrive at these quantities by various methods. In this pa- 
per, we attempt to examine in detail different theoretical mod- 
els and identify the sources of apparently large variations in 
galaxy merger rates between the models. We are interested 
both in defining the "en^or budget" of theoretical predictions, 
as well as in identifying the definitions, physics, and prescrip- 
tions that lead to substantial differences. 

In the follow sections we show that combining empirical 
halo occupation constraints with dark matter merger histories 
from high-resolution simulations leads to a relatively small 
(factor ^ 2) systematic uncertainty in merger rates (§|2]i. We 
then compare various different substructure based calcula- 
tions for subhalo-subhalo mergers and/or merger time delays 
(§|3]l. We highlight several caveats and potential problems in 
adopting these merger rates: without proper care, it is possible 
to obtain apparently order-of-magnitude different dark matter 
merger rates from the same simulation or analytic calculation. 

We reveal that most of the differences between model pre- 
dictions owe to the treatment of baryons in galaxy formation 
(§|4|i. Specifically, we show that cosmological hydrodynamic 
merger simulations lacking feedback do not reproduce the ob- 
served Mgai(Mhaio) relation (not just the normalization of this 
relation, but also its shape) and so do not map robustly be- 
tween halo-halo and galaxy-galaxy mergers (§ 14.1b . Semi- 
analytic models, on the other hand, perform well for cen- 
tral galaxies, but commonly have well-known problems re- 
producing observed properties of satellite galaxies; this ac- 
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counts for most of the discrepancies both between different 
semi-analytic models and between those models and empiri- 
cal approaches/calculations (§ 14. 2b . In §|5] we discuss the im- 
portance of consistent and appropriate mass ratio definitions, 
and the variations between different definitions, as well as the 
range in the effects of a given merger assuming changes in 
galaxy gas fractions, orbital parameters, and structural prop- 
erties. We summarize our conclusions and discuss future tests 
and improvements in § |6] 

Except where otherwise specified, we adopt a 
"concordance" cosmology with ( ^I m, ^a Ji, cr^, ns)= 
(0.3,0.7,0.7,0.9,1.0) and a IChabrierl (l2003h IMF, and 
appropriately normalize all observations and models shown. 
The choice of IMF systematically shifts the normalization 
of stellar masses herein, but does not otherwise change our 
comparisons. 

Throughout, we use the notation Mgai to denote the bary- 
onic (stellar-Hcold gas) mass of galaxies; the stellar, cold gas, 
and dark matter halo masses are denoted M,, Mgas, and Mhaio, 
respectively. When we refer to merger mass ratios, we use the 
same subscripts to denote the relevant masses used to define a 
mass ratio (e.g. /igai = Mgai.a/Mgai, i), always taken such that 

0<^l<l (Mgal.I >Mgal,2). 

2. OVERVIEW: COMPARING MODEL MERGER RATES 

The galaxy merger rate is fundamentally determined by two 
quantities: the halo/subhalo merger rate, and the manner in 
which galaxies populate halos (the halo occupation function). 
The halo mer ger r ate (the precise definition of which will be 
discussed in § 13. 2b is determined in cosmological simulations; 
there are two basic choices, however, for the halo occupation 
function. Semi-analytic models and cosmological hydrody- 
namic simulations attempt to predict this from first principles. 
Empirical models based on the halo occupation formalism 
simply adopt this quantity based on observational constraints. 

2.1. Methodologies 

Semi-analytic models and direct simulations are a well 
known and traditional means of modeling galaxy formation. 
In cosmological hydrodynamic simulations, galaxy formation 
is simulated from cosmological initial conditions. Cooling 
is computed self-consistently, but some prescriptions must be 
adopted for star formation (as a function of local gas prop- 
erties) and feedback from star formation and stellar evolu- 
tion as well as black hole (BH) and active galactic nuclei 
(AGN) fueling, growth and feedback. It is well-known how- 
ever that, in the absence of detailed prescriptions for both 
stellar and AGN feedback, such models do not reproduce the 
observed correlations between galaxy and halo mass - they 
overpredict the masses of both low and high mass galaxies, 
by a mass-dependent factor (i .e. predict a differen t magnitude 
and shape of Mgai(Mhaio); seelKeres et al . 2009a; Mailer et al.l 
120061; lOppenheimer et al.ll2010l;IChoi & Nagamin ell2009l) . As 
a consequence the most successful versions of such models 
have considered high-resolution "zoom-in" models of indi- 
vidual galaxies in formation, with extensive feedback pre- 
scriptions (e.g. jS ommer-Lars en et ani2003t iGovernato et al.l 
l2007HOkamoto et al. 2005; Sc annapieco et al.ll2008l) . Unfor- 
tunately, it remains prohibitively expensive to simulate cos- 
mological populations of galaxies (needed for the statistics to 
quantify galaxy-galaxy merger rates as a function of mass, 
mass ratio, redshift, and other properties) with the resolution 
and physical prescriptions necessary to follow these feedback 



and formation models, let alone to widely vary these pre- 
scriptions. And even many strong-feedbac k simulations over- 
predict the baryon conversion efficiency (Guo et al."2010at 
Abadi et al. 2003j_Okamoto et al. 2005; Governato et all 
2007; Scannapi eco et al.li2008t iPiontek & Ste inmetz 2009). 

Given this limitation, a more common approach to model- 
ing galaxy formation in an a priori fashion has been the con- 
struction of semi-analytic models. In these models, galaxy 
formation prescriptions are "painted onto" a dark matter back- 
ground. Dark matter halo locations, growth histories, merger 
rates, and other properties are taken from dark-matter only 
cosmological simulations and/or analytic extended Press- 
Schechter theory designed to match those simulations. Galax- 
ies are analytically assigned to each halo from the simulation, 
with their evolution followed in a Monte Carlo fashion ac- 
cording to a simple set of analytic prescriptions, designed to 
approximate the processes of cooling from halo gas, star for- 
mation, and stellar and AGN feedback. Such approaches have 
the advantage that they are inexpensive, allowing experimen- 
tation with a wide variety of formulations for e.g. feedback 
(with various free parameters), constrained so as to repro- 
duce observable quantities such as the galaxy mass function at 
z = 0. However, the prescriptions on which the results depend 
are approximations that usually must be applied to a variety of 
regimes well outside the range where they are observationally 
or physically (from e.g. direct simulations) calibrated. The re- 
sults in certain regimes can be unexpected or, in some cases, 
unphysical. In addition, it is well known that there are many 
degeneracies between small changes to different prescriptions 
- there is no unique solution to reproducing e.g. the z — 
stellar mass function, so the models are not guaranteed to be 
similar even if t hey are calibrated to the sam e observations 
initially (see e.g. lNeistein & WeinmannlbOlOl and references 
therein). 

In order to circumvent these physical uncertainties, 
models constructed using an empirical approach (the 
halo occupation formalism) have become popular as a 
means to predict the evolutionary paths of galaxies, in- 
cluding their clustering evolution, net growth, st a r for- 
mation histories, and merger rates (Zheng et al. '2007|; 
Yanetal. 2003; Conrov et al. 2007; Conrov & Wechsler 
2009; Perez-Gonzalez et al. 2008a; Tinker et al. 2005; Coora^ 
i2006: .Brown et al. 2008; Wetzel et al.. 2009: Stew art et al. 
l2009al;lHopkInse t al. 2010l: iMoster et alJl2010bl) . By match- 
ing the observed correlation functions of galaxies versus stel- 
lar mass, redshift, and other properties, one can empirically 
infer which halos the galaxies populate."' Independent ob- 
servational constraints on halo masses (discussed in detail in 
§ lU such as weak lensing, kinematic modeling. X-ray gas 
measurements, and group velocity dispersions yield comple- 
mentary constraints. Populating a cosmological simulation of 
halosH-subhalos in this manner, one can evolve it forward and 
determine where and when mergers occur. 

We should formally distinguish traditional strict HOD models from sub- 
halo abundance matching models. In the former, one paints galaxies onto host 
halos with no regard for halo merger histories or dark matter substructure, 
using clu stering as a constraint (see e.g. Yan et al. 2003; Yang et al. 200j; 
[Tinker et al.,2005; Kravtsov et al. 2004; Wech sler et al . 2006). In the latter, 
subhalo orbits are tracked in simulations and galaxies are assigned strictly 
to subhalos, but by matching the observed stellar mass function, so that the 
correlation function is a prediction of the model (which agrees very well; see 
e.g. Conroy et al. 2006; Vale & Ostriker 2006; Lee et al. 2009; Wetze l et alj 
|2()09H Wetzel & Whi tell20I0l ; lMoster et alj20I0bl ; lBehroozi et a l. 201(1). For 
our purposes however, the two methods give the same mass functions and 
clustering/substructure, so are treated interchangeably. 
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Such approaches are not, of course, models for galaxy for- 
mation - they populate galaxies in halos according to what 
is observed, and then (evolving them forward some arbitrarily 
small amount in time) can be used to predict quantities such as 
the merger rate (or what the merger rate must be if both dark- 
matter clustering simulations and the clustering observations 
are correct). They do not make any explicit statement about 
how the galaxies formed in the first place, nor the physics that 
are important. They require, for example, that the ratio of halo 
to stellar mass be high in low and high-mass systems, but in- 
clude no information about what feedback processes (believed 
to be stellar and AGN feedback, respectively) explain this. 
They are also, of course, limited by the observations: there are 
uncertainties in those observations corresponding to a range 
of Mgai(Mhaio) allowed, and there are redshifts and masses 
that the observ ations do not probe (see Mos ter et al. 2010^ 
iBehroozi et al.l i2010. and references therein), as well as de- 
generacies in fitting the model parameters that may hide crit- 
ical differences between different sub-populations (for exam- 
ple, those wh ich will or will not merge in a short time (see the 
discussion in lGuo_et al. 2010a). Roughly speaking, the ob- 
servational constraints themselves are robust above ^ 10''Mq 
at z = and 10" Mq at z ~ 2, but rapidly become more 
uncertain outside of this range. 

2.2. Semi-Empirical Models: Sources of Uncertainty 

With the advent of dark-matter simulations that can follow 
substructure at high resolution and tight observational con- 
straints on the stellar mass functions and clustering of galax- 
ies as a function of mass and redshift , semi-empi r ical ap- 
pro aches have become q uite robust. Stewart et al.l (l2009ah 
and lHopkins et alJ (1201 Ol) present predicted merger rates de- 
rived from such models, based on a variety of different as- 
sumptions and observational constraints. We refer to those 
pa pers for mode l details: the basic methodology is similar 
In lHopkins et al.l (12010.) . halo-halo merger hi stories are taken 
from the Millenium dark matter simulation dFakhouri & Mai 
20081) (these are updated with the Millenium-II results in 
Fakhou ri et al. 2010, but the results for our purposes are iden- 
tical), and in S tewart et^ (|2009a) from independent dark 
matter simulations "; each halo is then populated with galax- 
ies using the abundance-matching methodology (although 
variations and alternative halo occupation approaches are at- 
tempted), ensuring (by construction) a match to observed 
galaxy mass functions and correlation functions at all red- 
shifts. In short time intervals, the model halos are evolved 
forward; when halos merge, either the subhalos are followed 
until they are no longer resolved or a "dynamical friction 
timescal e" is assigned as calibra ted from high-resolution sim- 
ulations (Boylan-Kolchin et al. 2008), and at the end of this 
time the hosted galaxies are assumed to merge. Summing over 
the entire simulation volume, this gives the predicted galaxy- 
galaxy merger rate. 

The two independent models find predicted merger rates 
that agree reasonably well: typical variations in the model in- 
puts lead to factor < 2 changes in the predicted merger rates, 
at least at redshifts z < 2 where the halo occupation function is 
reasonably well-constrained by observations. From this per- 



spective, the advantage of the semi-empirical models is that 
they do not attempt to model galaxy formation - as a result, 
differences between various predictions must owe either to the 
halo occupation distribution adopted (i.e. constraints on the 
galaxy mass-halo mass distribution), or to the "background" 
dark matter dynamics. 

Figure [T] summarizes the comparisons from iHopkins et al] 
(i2010i) .'^ The different model choices are discussed therein, 
and we will break each down in detail in what follows. 
For now, we simply illustrate the different effects on the 
merger rate.'^ Specifically, we compare the predicted ma- 
jor merger rate (mergers of mass ratio /igai > 1/3, where 
^gai = Msecondary/A^primaiy < 1 in tcrms of the galaxy baryonic 
mass Mgai) for intermediate-mass galax ies as a function o f 
redshift. We begin with the model from Hopk ins et al.l (1201 Ol) . 
and consider the effects of the choices below. We take this 
as our basic HOD or semi-empirical m odel throughout, fre- 
quently comparing with the model from iStewart et alj (l2009ah 
(which uses different sets of empirical HOD constraints, dark 
matter simulations, and methods of subhalo tracking to iden- 
tify mergers), but with many of the variations below consid- 
ered separately (each can be considere d its own mode l). 

(1) Cosmology: We re-run the IHopkins et aP (1201 Ol) 
model adopting the c osmological parameters from WM APl 
dSpergel et al. 2003), WMAP3 (ISoergel et alJ l2007l) and 
WMAP5 ( Komatsu et al. 2009), as well as a "concordance" 
model with (f^M, f^A, /i, erg, n.)=(0.3, 0.7, 0.7, 0.9, 1.0). Be- 
cause the halo occupation statistics are constrained to repro- 
duce the same stellar mass function, most of the difference 
between cosmologies (for example, in the predicted halo mass 
function) is effectively normalized out in this approach. 

(2) Halo-Halo Merger Rates: The iHopkins et alj 
(1201 Ol) model a dopts the halo-halo merger rate fits from 
iFakhouri et al.l (I201C ), from the Millenium simulation 
('Sprin gel et al.l l2b05b ). An alternative dark-matter simu- 
lation, of comparable resolution, with halo merger rates 

determined u s ing a d ifferent methodolog is described in 
'Stewa rt et all ( l2009al). Another is fou nd in Gottlober et ajj 
(2001) (see also iKravtsov et all 120041: fZentner et al. 2005); 
they quantify the fit separately to field, group, and cluster 
environments. We also consider the galaxy-galaxy merger 
rates directly determined fr o m cos mological hydrodynamic 
simulations in ' Mailer et alJ ([200^. It is well-known that 
such simulations yield galaxies that are overmassive relative 
to those observed (i.e. predict an HOD in conflict with that 
observed) - therefore we only use these mergers as "markers" 
of where the galaxy-galaxy mergers should occur, and (in 
order to study the stellar mass and mass ratio dependencies) 
re-populate the galaxies with the "correct" masses according 
to the HOD. 

(3) Substructure: Instead of using a halo-halo merger 
rate with some "delay" applied to model subhalos in the 
HOD, we can attempt to follow subhalos directly after the 
halo-halo merger, and define the galaxy-galaxy merger when 
the subhalos a re fully merged/destroyed. We compare the 
IHopkins et al.l ([2010) rates to those obtained tracking the 
halo-H subhalo popula tions in cosmological simulations from 
IStewart et al.l (l2009al) (populating subhalos according to the 



' ' The simulation is run with the ART code f Kravt sov et alj[T997l) . with 
distinct halo identification and subhalo tracking algorithms from those used 
in the Millenium simulations iKravtsov et al. 2004), and a different subhalo 
'destruction' criteria based on the fractional mass loss via stripping (subhalos 
merged when this falls below a fixed fraction of their inflow mass) as opposed 
to an absolute mass resolution limit. 



The merger rates derived from this model, including several of the 
variations discussed here, can be obtained from the "merger rate calculator" 
routine publicly available at http : / /www . cf a . harvard. eciu/~phopklns/ 
Site/mergercalc .html. 

Henceforth, we take the term "merger rate" to refer to the galaxy-galaxy 
merger rate, unless otherwise specified. 
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Fig. 1 . — Comparison of the galaxy-galaxy major merger rate (mergers per galaxy per Gyr) for Mg-^i ~ lO'" — lO" Mq galaxies as a function of redshift, 
from different variations in semi-empirical models - i.e. mod els in which g alaxies are placed in halos strictly according to observational constraints. Black line 
in each case corresponds to the "default" model from Hopki ns et al.l )201(B . Red dotted range is the range allowed by observations compiled in that paper Top 
Left: Changing the halo occupation model (observational constraints used to place galaxies in halos); each of the model choices listed is discussed in detail in the 
text (§ 12.2^ The choices shown bracket the range allowed by a number of independent observational constraints. Top Center: Changing the merger "timescale" 
- i.e. timescale used to estimate the delay between a halo-halo merger and the resulting galaxy-galaxy merger (during which the satellite no longer grows). Top 
Right: Directly following subhalos in simulations, and assigning mergers when those subhalos are destroyed/fully merge into the central halo, instead of using 
the merger "timescale" approach. Results are shown from various methods of following subhalos in simulations and EPS trees. Bottom Left: Adopting estimates 
of the parent halo-halo merger rate, from different dark matter simulations. For the merger rates from Gottlober et al. (2001), we also show the results using their 
merger rates in different environments (labeled). We also compare rates using cosmological hydrodynamic simulations to "tag" the halo-halo mergers of interest 
by identifying their resulting galaxy-galaxy mergers. Bottom Center: Changing the cosmological parameters. The allowed variations in these semi-empirical or 
HOD models lead to factor ~ 2 — 3 differences in the predicted merger rate. 



same HOD). We also compare with the results using the dif- 
ferent subhalo-based methodology described in'Hopki ns et al.l 
(|2008b) (essentially, beginning from the subhalo mass func- 
tion constructed from cosmological simulations and evolv- 
ing this forward in short time intervals after populating it, 
at each time, according to the HOD constraints). We com- 
pare two different constructions of the subhalo mass func- 
tions: that frorn cosm ological d ark-matter only sim ulations in 
iKravtsov et"!!] (2004) (see also IZentner et al.ll200 5) and that 
from the extended Press-Schechter formalism coupled to ba- 
si c prescriptions for subhalo dynamical evolution, described 
in Ivan den Bosch et all (l2005h . 

(4) Merger Timescales: Another commonly applied method 
in determining galaxy-galaxy merger rates as distinct from 
halo-halo mergers is to assign a simple delay to the galaxy- 
galaxy merger after the halo-halo merger, based on e.g. the dy- 
namical friction time for the subhalo-i-galaxy orbit to decay, as 
calibrated in high-resolution N-body simulations that resolve 
the galaxy a nd follow the entire system self-consistently. The 
iHo pkin s et al. (2 010.) model adopts the traditional dynami- 
cal friction time, using the calib r ation f rom numerical sim- 
ulations in iBovlan-Kolchin et al.l (l2008l) . We compare with 



the model re-run using the characteristic timescale for pair- 
pair gravitational capture in group environments, calibrated to 
simulations following Mamon (2006). We again repeat, con- 
sidering capture in angula r momentum spa c e rathe r than grav- 
itationally, following Binnev & Tremaiii^ (Il987l) . as well as 
the gravitational cross-section timescale (similar to the group 
capture timescale) for capture between passages in e.g. loose 
group or field environments, ca librated from simulations in 
iKrivitsky & Kontorovichid 19971) . 

(5) Halo Occupation: We re-run the IHopkins et al.l(l2010l) 
model using a different set of halo occupation constraints 
to determi ne the ga l axy m ass-halo mass relations. The de- 
fault Hopk ins et al.l (|2010|) employ s the fi ts to these dis- 
tributions from I Conroy & Wechsl^ (l2009l) . which adopt a 
monotonic rank-ordering approach to assign galaxies to ha- 
los in a manner constructed to match observed mass func- 
tions and clustering ( see the cornpilation in Fig. 1 therein 
and Bell et al. 2003 bl iDrory et all 120051: iBorch et all 120061: 
Fonta na et al.li2006c IPanter et al.ll2007l: iPerez-Gonzalez et alj 
2008b). Gas masses are added to galaxies according to the 
observed correlation between galaxy gas frac tions and stellar 
mass, from observations spanning z = 0-3 dBell & de JongI 
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2000HMcGaughl2005HCalura et al.l2008l:[Shaplev et al.l2005 



Erb et alJ 120061: iPuech et alJ 120081: Mannucci et all 1 2009 



Cresci et al.' '2009^, ' Forster Schreiber et al.' '2009^, 'Erb* '2008. 
Mannucci et al. 200^ compiled in popkins et al. (2010) 
(§ 2.2.2 therein). We compare to the results using the same 
monotonic rank-ordering HOD method ology, but with galaxy 
stellar m as s function constraint s from iPerez-Gonzalez et alj 
(l2008bl) or lFontana et al.l (l2006 l) (each determined up to red- 
shift z^4). We also compare to the results selecting the fit to 
M, (Mhaio) and its scatter for central and satellite galaxies from 
the observed SDSS clustering at z = (Wan g et al. 2006 ): in 
this case, we simply adopt the z = fit at all redshifts - we do 
not allow for evolution. 

We will discuss each of these aspects in detail in what fol- 
lows. However, it should already be clear in Figure [T] that 
these choices, while not negligible, introduce only factor ^ 2 
uncertainties in merger rates. Order-of-magnitude effects, it 
seems, must depend either on some particular combination 
of effects above, or relate to the baryonic physics of galaxy 
formation, such that the galaxy-halo mappings predicted are 
significantly different from those empirically determined by 
HOD models. 

2.3. A Priori Models Compared 

We now consider various a priori models for galaxy for- 
mation. In these cases, there are of course similar choices 
to those discussed above, for quantities such as the halo-halo 
merger rate, merger timescales and/or subhalo tracking, and 
cosmological parameters. However, the halo occupation dis- 
tribution is not determined from observations, but is predicted 
as a consequence of attempting to model galaxy formation in 
an a priori manner based on some set of physical assumptions. 
That is not to say this is completely unconstrained, however: 
the models (especially semi-analytic models) are generally 
adjusted so as to reproduce various observations such as the 
stellar mass function and large-scale galaxy clustering. As a 
result, at least certain portions of the model should be similar 
to the HOD models. 

The specific models we will consider throughout the paper 
are briefly as follows: 



(a) ISomerville et alj dlOOSal): A semi - analyti c 
model updating those in ISomerville & Primackl (Il999l) : 
ISomerville et al.l (120011) . Halo growth is tracked using the 
extended Press-Schechter formalism (using a slightly modi- 
fied version of the method of Somerville & Kolatt ( 1999), as 
described in Somerv ille et al.l (l2008al) . to match the results of 
A^-body simulations), with merging of sub-halos in virialized 
halos followed with appropriate dynamical friction times 
(including mass loss and tidal destruction). Cooling is cal- 
culated in the typical manner, distinguishing between "cold" 
and "hot" mode regimes (gas accretes onto central galaxies 
in a dynamical time in halos with Mhaio ^ IO'^Mq, but forms 
quasi-static pressure-supported halos at larger masses, as 
motivated by cosmological simulations; see e.g. 'Ker es et aLl 
12005). Disks form assuming conservation of s pecific angular 
momentum as describe d in IS omervill e et al.l (|2008b), and 
star formation follows iKenni cutt (19981) with supernova 
feedback. Black holes and bulges grow in mergers, with 
feedback associated with bright quasars and low-luminosity 
AGN, the latter of which suppresses new cooling in the "hot 
mode" of massive halos. This model adopts t he prescription 
for bulg e formation calibrated to simulations in lHopkins et al.l 
( l2009al) . given the gas fractions and mass ratios of mergers. 



Satellite galaxies are assumed to lose their hot gas reservoirs 
to their new parent halos (and cease to accrete new gas) upon 
halo-halo merger in the EPS tree, but do not lose the recycled 
gas from star formation internally in their pre-existing cold 
disks. 



(b) [Khochfaii&SnH (l2009l) (see also WeinzirLeLafl 
2009h: This s emi-ana lytic model differs from that of 



Somerville et al.l (l2008al) primarily in that there is no "quasar 



mode" of feedback, and the prescriptions for the low- 
luminosity "radio mode" of AGN feedback and stellar wind 
feedback differ somewhat in detail. Galaxy structural prop- 
erties are followed using e nergetic argu ments and scalings 
motivated by simulations dKhochfar & Silk 2006). Bulge 
formation is followed with a more simplified recipe common 
in most semi-analytic models: the secondary stellar disk is 
destroyed (added to the bulge), but the primary disk is only 
violently relaxed in a m ajor merge r , and s tarbursts in mergers 
are tracked following ICox et all ( 120081) . The differences 
stemming from the bulge formation model are discussed in 
detail in Hopkins et al. (2009b). Satellite galaxies retain their 
cold and hot gas reservoirs until the merger with the primary 
central galaxy (but do not accrete new gas). 



(c) iBower et aP (120061) : This iteration o f the Durham 
Galform SAM revises previous versions (iLacey & Cold 



Il993t iBenson et al.1 120031 iBaugh et al.l l2005h by adopting 
halo merger trees determined in N-body simulations, sp ecifi- 
cally the Millenium simulation (ISpringel et al.l l2005bl) . and 
tracking substructure within those halos, as well as adding 
a mode of low-luminosity AGN feedback (a "radio mode") 
which suppresses new cooling in "hot mode" halos. Cooling 
is followed as in the models above, and bulge format i on with 
the prescription later adopted by iKhochfar & Silkl (l2009l) . 
with the addition of allowed bulge formation in a strong 
disk instability mode that dominates at high redshift. This 
mode sets in when the disk becomes self-gravitating, and 
in consequence the disk is assumed to largely collapse to a 
bulge. Satellite galaxies are assumed to lose all of their hot 
gas to the new parent halo upon first linkage of the parent 
halos (often several /?vii), and gas which galaxies "expel" 
owing to the assumed stellar wind mode of feedback is 
considered part of the halo gas and treated similarly. 



(d) iFont et alj (|2008^: The most recent ite ration of the 
Durham SAM . The model improves on that of iBower et al.l 
(l2006l) by including a physically-motivated prescription for 
the removal of hot gas from subhalos after they are first 
linked to the main halo group (as a function of e.g. their 
position in the primary halo and the background gas density; 
motivated by ram-pressure stripping models). It is otherwise 
identical. We will show in § l4.2l that this has a large effect on 
the predicted merger rates. 



(e) ICroton et alJ (120061) : Another semi-analytic model 
based on the out put of the Mi l leniu m simulation and broadly 
similar to that in Bower et al. (2006), but with important dif- 
ferences in the prescriptions for stellar feedback and cooling 
in both central and satellite galaxies (see § 14.2b . Accretion 
in the "cold mode" is delayed by a gal axy dynamic al time. 
Satellite halos are instantly stripped as in lBower et alJ (12006): 
the mass expelled by star formation is contained in a separate 
reservoir to be recycled, but this too is lost to the parent halo 
for satelUte galaxies. The merger trees are constructed from 
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the same simulation, but with different algorithms and cuts. 
The model also includes a much weaker disk instabili t y mod e 
of bulge formation, relative to that in iBower et alj (l2006l) . 
The implementation of AGN feedback is quaUta tively similar 
between the two models, but quantitatively the Croto n et alJ 
(12006.) mode is "stronger" and yields more radio power. 



(f) Ide Lucia & BlaizotI dlOOTl) : The most recent iteration 
of t he Munich SAM, th is is a somewhat revised version of 
the ICroton et al.1 (|2006|) SAM. The primary modifications 
include the adoption of a different stellar IMF and revised 
dust attenuation calculation, and the removal of a dynamical 
time delay in "cold mode" accretion (making "quenching" 
more rapid). Also, in this version, dynamical friction times 
assigned to galaxies before mergers are increased relative to 
the canonical dynamical friction time by a factor ^ 2, with 
the goal of suppressing mergers at the massive end of the 
galaxy mass function (since the models tend to have difficulty 
reproducing this). 



(g) lMaUer et al J (12006): The results of direct cosmological 
hydrodynamic simulations, incorporating cooling and star 
formation self-consistently, but without feedback of any 
kind from stars or BHs. Galaxy-galaxy mergers are directly 
identified, and their mass ratio quantified at the time just 
before merger Hydrodynamic simulations with very large 
dynamic ranges in mass and redshift remain prohibitively 
expensive, so the comparison is restricted to intermediate and 
high-mass galaxies at redshifts z < 2. As the authors note, 
without feedback, the mass function at z = is over-produced 
(and the break is not as strong as observed). They adopt a 
mean correction for this; essentially, re-normalizing their 
stellar masses by a systematic factor of ^ 3. 



Figure |2] compares the predicted merger rates from these 
models. We compare these with observational determinations 
of the galaxy-galaxy merger rate, from close pairs. Specif- 
ically, from the fraction of major (/igai > 1/3) parrs with 
small projected separations r^, < 20/z^'kpc (often with the 
additional requirement of a line-of-sight velocity separation 
< 500 kms^'), and stellar masses in a specific range (labeled). 
For each mass bin, the pair fractions as a function of redshift 
can be empirically conver ted to a merger rate using the merger 
timescales at each radius. iLotz et alj ( l200 8a) specifically cal- 
ibrate these timescales for the same projected separation and 
velocity selection from a detailed study of a large suite of hy- 
drodynamic merger simulations (including a range of galaxy 
masses, orbital parameters, gas fractions and star formation 
rates) using mock images obtained by applying realistic ra- 
diative transfer models, with the identical observational crite- 
ria to classify mock observations of the galaxies at all times 
and sightlines during their evolution. For this specific pair 
selection criterion, they find a median merger timescale of 
Wigei ~ 0.35 Gyr, with relatively small scatter and very little 
dependence on simulation parameters (±0.15 Gyr). We use 

'* The merger timescale from simulations at this radius is somewhat 
shorter than the time obtained assuming dynamical fiiction and circular orbits 
in e.g. an isothermal sphere, as has commonly been done (this is assumed in 
e.g. both Patton et al. 2002; Kitzbichler & White 2008). This owes to two ef- 
fects: first, angular momentum loss at these radii is not dominated by dynam- 
ical friction (at least in the traditional sense of a small mass moving through 
a smooth, isotropic velocity dispersion background), but rather by exchange 
in strong resonances between the baryonic components that act much more 



this median fmerger to convert the observations to a merger rate. 
Completeness corrections are discussed in th e various papers; 
we also adopt the standard correction from Patton & Atfieldl 
(1^08), calibrated to high-resolution simulations, for the frac- 
tion of systems on early or non-merging passages (to prevent 
double-counting systems on multipl e passages); but t his is 
relatively small (20 - 40%; see also iLotz et al.ll2008ah . For 
more detailed discussion of these observations and related 
constraints from othe r methods and observations, we refer to 
iHopkins et all ( 120101) . 

We are not attempting to comprehensively discuss the pre- 
dicted merger rates here: those rates are presented in other 
papers in much greater detail.'^ However, we do wish to ex- 
amine what can affect the predicted merger rate. It is clear in 
Figure |2] that the variation in predictions between different a 
priori models is significantly larger than that within the semi- 
empirical models. This motivates our further examination of 
the differences and similarities in these models. 

For clarity, we have limited the number of semi-analytic 
models shown in Figure |2] However, we have compared with 
a number o f other s uch models and find similar results (in- 
cluding e.g.lBenson et al. 2003; Baugh et al. 2005; Kang et aj] 
20051; ICattaneo et al.ll2006; Monaco et al.,,2007; .Bertone et all 



20071). The models considered span the representative behav- 
ior in these other semi-analytic models as well. 

3. THE DARK SIDE 

First, we examine the differences in the dark matter side of 
the merger rate, i.e. effects on the predicted merger rates that 
are independent of the detailed baryonic physics of galaxy for- 
mation. 

3.1. Cosmological Parameters 

As shown in Figure [1] the assumed cosmology makes lit- 
tle difference for merger rates, so long as models are nor- 
malized to produce similar stellar mass functions at each red- 
shift. Present constraints on the cosmological parameters are 
also quite strong (Komatsu et al. 2009); there is little freedom 
to change the cosmology in a significant fashion. At high 
redshifts, the differences as a consequence of e.g. the value 
of (Ts become larger, but so do the other sources of uncer- 
tainty below - at no point do we find that the uncertainties in 
cosmological parameters dominate the variations in predicted 
merger rates. Moreover, most of the models considered here 
adopt an identical cosmology (for example, all of the semi- 
analytic models in Figure [T] have been run or have available 
versions adopting the "concordance" cosmological parame- 
ters defined in § [TJ, and exhibit nearly identical variation be- 
tween models. 

efficiently. Second, by these radii, even initially circular orbits have become 
highly radial, leading to shorter merger times. Because of these effects, the 
remaining merger time at this scale depends only weakly on initial condi- 
tions or orbital parameters - essentially, these processes have erased most of 
the "memory" of the original orbital configuration. 

For predicted merger rates from the Somervl ile et alj (2008aD mo del, 
see Hopkins etal. (2009b, 2013). For those from Khochfar & STffl 
(2009), see Weinzirl et al. (2005?^^ From power et al. (2006), see MateiJ 
(2008) and PaiTy et al. (2009). From de Luci a & BlaizotI <2007ft. see 
Guo& White (2008), Kitzbichler & White (20081), and IBertone et"an 
^2007). The Mailer et al. (2006) rates are discussed therein. The model 
predictions from de Lucia & Blaizot (2007); Bower et al. (2006); Font et al] 
12003) can also be obtained online from the Millenium public database, 

at [http : / /www . mpa-garching . mpg . de/millenn Tum/l or 

http://icc. dur .ac.uk/ index . php?content=MiIlenn ium/MillenniuiTi[ qi 
shown here from these two models are obtained from the same database). 
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Fig. 2. — As Figure[T] but comparing the HOD-based models with different predictions from a priori galaxy formation models, in different galaxy stellar 
mass interv als ( log (Mn /Mp^), as labeled). We compare the standard semi-empirical HOD models from Hopkins et al. (2010) and Stewart et al. l2009a) (de- 
scribed in § 12. 2) to the direct predictions of cosmological hydrodynamic simulations (with cooling and star formation, but without feedback from either stars or 
AGN; Mailer et al. 2006) and the predictions of various semi-analytic galaxy for matio n models (Croton et al. 2006; Bower et al. 2006; de Lucia & Blaizot 200'^; 
Ehochfar & Silk 2009; So merville et al. .2008a ; Font et al. 2008, described in §|23). We compare these predictions with observations in each mass interval, 
based on close pair counts (r < 20A~' kpc) with calibration of merger timescales specific for each observed sepa r ation and selection crite ria from lar ge suites 
of high-resolution galaxy merger simulations ^Lotz et al. 2008a). The observations are compiled in Hop kins et alj j2010l) from Kartaltepe et alj j2007l blue tri- 
angles), Conselice et al. (2009, blue inverted triangles), Lin et al. (2004, 2008, pink circles), Xu et al. ( 20041 blue circle). IDe Propris et alj j2005l black asterisk), 
iBluck et aL 42009. blue squares), .Bundy et al., t2009» green stars), and iBell et aljt2006bia. red pentagons). 



Comparing the halo mass function and halo accretion 
histories with diffe rent cosmological parameters (see e.g. 
iNeistein et al.ll2006l) . the dominant effect is the predicted halo 
mass function shifting to higher masses with larger erg. How- 
ever, if we re-normalize the model to enforce the same ob- 
served galaxy stellar and baryonic mass functions and clus- 
tering, then these differences are largely normalized out, at 
least for variation in erg in the relatively small observationally 
allowed range (here 0.7 — 0.9 considered; much larger differ- 
ences may significantly change dynamical friction times and 
other higher-order effects). Elahi et al. (2009) show that the 
quantity of greatest importance for our conclusions, the nor- 
malized substructure mass function or (equivalently) dimen- 
sionless merger rate (mergers per halo per Hubble time per 
unit mass ratio) is almost completely independent of cosmo- 
logical parameters including e.g. the power spectrum shape 
«s ^ 1 and amplitude over the range of variations here (not 
until one goes to much larger effective ^ 3 does one see 
this function change shape). 

3.2. Halo-Halo Merger Rates 

3.2.L Comparison of Different Simulations and EPS 
Calculations 

Even given the same cosmological parameters, different 
methods used to determine the halo-halo merger rate can 
yield different answers. In Figure [3] we plot the depen- 
dence of the halo-halo merger rate on halo mass ratio /ihaio = 
Mh-dio.i/Miinio.i < 1 and redshift, from different cosmological 
simulations and analytic calculations. Specifically, we plot 
the number of mergers per unit redshift above some /ihaio at 
z = 0, and the evolution with redshift for major (/ihaio > 1/3) 
mergers (evaluated in a narrow mass interval around lO'^M© 
halos, in terms of post-merger mass, although all these models 
find that in these units, the mass dependence is weak). 

We consider s everal different determ inations of the halo- 
halo merger rate. iFakhouri & Mai ( l2008l) calculate the merger 



rate in the Millenium simulation in the usual fashion (deter- 
mining the number of mergers per descendant halo in terms of 
the halo friends-of-friends group mass and mass ratio /ihaio < 
1). These are updated with the results of th e higher-resolution 
Millenium-II simulation in Fakhouri et al. (2010), but the re- 
sults are for our purposes identical. Genel et al. (2009) ana- 
lyze the same simulation, but determine the merger rate us- 
ing another methodology: defining halo masses in a differ- 
ent manner and quoting the merger rate per progenitor halo; 
they discuss in detail how this leads to subtle but potentially 
important variations in the merger rate fits. However, con- 
verted to the same units (following the methodology in their 
Appendix), the results are similar Stewart et alj (12009 a l) de- 
termine the halo-halo merger rate from an independent dark 
matter simulation, using a slightly revised methodology and 
different time-step criteria. We also compare these simulation 
results with an analytic app roximation: using the methodol- 
ogy in Ivan den Bosch et aU (l2005i) to predict the merger rate 
analytically from a modified extended Press-Schechter for- 
malism (note that these authors consider only the main branch 
mergers; because we are considering instantaneous merger 
rates here, this makes no difference, but it will lead to a dif- 
ferent integrated numbe r of mergers if not accounted for). 
INeistein & Dekell (l2008l) also present Press-Schechter trees 
modified to better match numerical simulations. Yet another 
methodology is described in Hopkins et al. (2008b); those au- 
thors take the subhalo mass function from cosmological dark 
matter simulations, and infer the halo-halo merger rate re- 
quired to maintain this mass function allowing for dynami- 
cal evolution of the subhalos (i.e. assume steady-state subhalo 
mass functions over a narrow redshift interval, given certain 
initial conditions). We specifically adopt th eir estimate us- 
ing th e su bhalo mass functions presented in iKravtsov et alj 
(12004 and lZentner et al.1 (l2005l) . In that paper, however, they 
consider several other simulation results as well, and find sim- 
ilar results in all cas es (Shaw et al. 2006; Sprinael et al. 200J; 
iTormen et all l2004t lOe Lucia et al.1 12004t iGao et al.l 12004 : 
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Fig. 3. — Left: DM halo-halo merger rates (mergers per halo per unit redshift) as a function of halo mass ratio (at z = 0), from different simulations 
and methodologies shown. Fakhouri & Ma (2008) and Genel et al. (2009) extract merger rates from the Milleniu m and Millenium- II DM-only simulations, 
with different methodologies. Stewart et al. 1 2009^) determine rates from independent DM simulations. The Ki'avtso v et alj {2004) rates are based on the 
integral constraint from subhalo mass functions, rather than an instantaneous merg er r ate. R ates from Mailer et al. 1 2006) are from hydrodynamic cosmological 
simulations (using galaxies as halo tracers). Rates from lvan den Bosch et alj 120051) and lNeistein & Dekel (^2008) are constructed with the EPS formalism. Right: 
Corresponding halo major merger rate versus redshift. Most different simulated merger rates agree well (within ~ 20 — 50% for > 0.1), provided definitions 
are appropriately accounted for. Results are shown for M^aio = IO'^Mq (at the given redshift), but depend very weakly on halo mass. 



iNurmi et al.lf2006l) . 

The mergers of halo barycenters are the real mergers of in- 
terest, and it is possible that the presence of baryons could 
alter halo mass functions and merger rates at some level. 
These and other baryonic effects can only be followed in cos- 
mological hydrodynamic simulations. Of course, such mod- 
els may not reproduce the correct observed galaxy baryonic 
masses. However, we can still use such simulations as trac- 
ers of galaxy-galaxy mergers; using the information about the 
position and mergers of galaxies, but replacing the galaxy 
masses according to the observed HOD (essentially using 
galaxies in the simulation to 'tag' the relevant halos merg- 
ing).'^ This also allows for the possibility that the baryonic 
material changes the halo-halo merger rate, by e.g. contribut- 
ing to halo self-gravity, although this is typically estimated to 
be a ^ 10 — 20% effect for the m ajor-merger regime of inter- 
est (see e.g. [Weinberg et all2008l) . Figure [3] shows the results 
of this procedure given the cosmological smoothed particle 
hydrodynamics simulations in Mailer et al. (2006). 

Overall, there is quite good agreement between these dif- 
ferent calculations - halo-halo merger rates as determined in 
cosmological simulations are generally converged, at least for 
intermediate masses and redshifts z < 2 — 3, at the ^ 50% 
level or better. The scatter in different simulation-based de- 
terminations grows to a factor ~ 2 at high/low masses. 

Furthermore, many of these differences owe to choices of 
definition that, treated properly, w ill not enter i nto the galaxy- 
galaxy merger rate. For example. iGenel et"an (l2009h demon- 
strate in detail how changing the halo mass definition and def- 
inition of merger rate in terms of mergers per progenitor (and 
allowing for /ihaio > 1; i.e. mergers into larger halos) leads 

We do this by simply replacing the rank-ordered galaxy masses with 
those obtained by abundance matching. 



to differences in the halo-hal o merg er rate relative to the fits 
presented in iFakhouri & Mai (120081) & iFakhouri et al] (1201 Oh 
from the same simulation, and the differences can be signif- 
icant (factor ^ 2). However, the merger trees in terms of 
e.g. halo barycenters in both cases are almost identical. In 
terms of galaxy-galaxy merger rates, the precise halo mass 
definition, for example, should make no difference, so long as 
the halo occupation statistics are appropriately renormalized 
to yield the same galaxy mass function and galaxy cluster- 
ing versus luminosity (halo mass is simply a label to use in 
tagging galaxies to halos). Many of the differences in Fig- 
ure [3] therefore, are actually reduced in terms of their effect 
on galaxy-galaxy merger rates, so long as sufficient care is 
taken in correcting for the definitions used. 

Considerable effort has also gone into tuning EPS -predicted 
merger rates to bett er correspond to the results of N-body 
simulations (see e.g. IZentne^l2007l : iNeistein & De kel 2008). 
With these improvements, many current-generation EPS trees 
differ from the simulation results they are tuned to by as little 
as ^ 10 — 20%. For detailed c omparison o f EPS and simu- 
lation merger rates, we refer to iGenel et al.l (12009). Clearly, 
this is not the dominant source of uncertainty in merger rates, 
as the differences between various counts of mergers from the 
same simulation differ by larger amounts, depending on their 
exact methodology. 

3.2.2. Some Caveats in Adopting Halo-Halo Merger Rates 

We do not mean to imply that choices of definition are en- 
tirely negligible. Depending on the application, a halo merger 
rate defined in terms of e.g. progenitor galaxies or descendant 
galaxies is more appropriate (the former being an estimate of 
the probability that a given halo will have a merger; the latter 
being an estimate of the probability that a given halo did have 
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Fig. 4. — As Figure|3] showing the differential number of halo-halo merg- 
ers per lo g mass ratio (at z = 0) from the Millenium DM-only simulation. 
The lFakhouri & Ma qOO S ) and Genel et al., (,2009) rates are intended to rep- 
resent the halo-halo merger rate. We compare with the instantan eous rate 
taken directly from the Durham semi-analytic model halo trees I Marker et al.l 
12006). or from a simple query of the simulation database of the Munich halo 
trees using the instantaneous mass ratios just before the secondary is "losf 
in the tree. Despite this being the same set of trees, these "instantaneous" 
mass ratio definitions lead to an apparent large suppression of major halo- 
halo mergers. This is an artificial effect of the definitions used; we show the 
differences to illustrate the importance of accounting for different halo mass 
ratio definitions in constructing the halo-merger rate. 
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Fig. 5. — As Figure [5] comparing the halo merger rate as a function of 
redshift from the Fakhouri & Ma 1 2008) trees. We compare the merger rate 
using the minimum available time spacing ( At = 0.00 Gyr, solid; in fact the 
spacing is finite but the halos are tracked on-the-fly in the simulation, so it is 
small for our purposes), to the rate inferred from comparing snapshots only 
with the given time spacing, if the mass ratio is defined by comparison of 
subsequent snapshots. Large timestepping means major mergers appear mi- 
nor or are lost. If a quantity depends super-linearly on halo mass (as e.g. 
galaxy mass does at low masses), then the apparent suppression will be en- 
hanced. 



a merger). And there are means of constructing halo catalogs 
that could give quite different halo-halo merger rates, despite 
being ultimately a reflection of the same merger tree. The ap- 
plication of any halo merger rate requires careful accounting 
for consequences of such definitions. 

Figure |4] illustrates this. We compare the halo-halo 
merger rates constructed fro m the Milleniu m simula- 
tion by Eak houri & Ma (20081; iFakhouri et aP (2010) and 
iGenel et aL (2009) to those constructed from the s ame simula- 
tion for use in the Durh am semi-analytic models tBower et al.l 
I2OO6I: iFont et al1l2008l) . according to the method described in 



iHarker et all (l2006h . There is an abrupt cutoff in the number 
of major halo-halo mergers in the latter tree. The differences 
are an alyzed in detail in the Appendix of IFakhouri & Mai 
( I2OO8I) . but they largely owe to two post-processing cuts made 
after the original halo-i-subhalo tree construction: subhalos are 
removed from their parent friends-of-friends (FOF) identified 
halo if their centers are outside twice the half-mass radius, or 
if they have retained at least 75% of their mass at the time they 
were last an independent halo. Subhalos (and their progeni- 
tors) that are removed from groups in this way are discarded 
from the trees. 

The intention of these cuts is to reduce the number of spuri- 
ous linkings from the FOF group finder, and both cuts above 
probably do so. However, they also suppress the number of 
halo-halo major mergers in the tree for three reasons. First, 
the discarded subhalos change the shape of the halo mass 
function - some high-mass halos lose a non-negligible frac- 
tion of their mass and, on the exponential tail of the mass func- 
tion, this can have large effects (since small changes in mass 
correspond to large changes in number densities). Second, the 
secondary (sub)halo in a major merger is (fractionally) more 
robust to stripping than a small secondary in a minor merger, 
so major mergers are preferentially affected by the mass re- 
tention threshold. And third, in an equal mass merger, by 
definition the halos first overlap when the secondary center is 
at ~ 2Ryt from the primary (in fact given the FOF algorithm 
used for linking, the distance can be 3 —5Rvt), so many 
of its subhalos will be instantaneously split off (being outside 
twice the half-mass radius at this instant) and thus not con- 
tribute to the secondary halo mass at the moment of merger 
(so the instantaneous mass ratio suddenly "dips" just before 
the merger). 

This is not necessarily problematic for galaxy merger rates. 
The trees so constructed are designed for tracking the galaxy 
population in a specific semi-analytic model; they are not de- 
signed to represent the global halo-halo merger rate. The ef- 
fects above mainly pertain to definitions; they mostly apply to 
the early stages of subhalo-halo mergers, not necessarily the 
final, relaxed stages when a galaxy merger will occur. As such 
they only indirectly affect galaxies, for example, via the small 
changes to the total halo mass introduced and corresponding 
gas mass assumed available for cooling. So the apparent dif- 
ference in halo-halo merger rates translates much less directly 
in terms of galaxy-galaxy merger rates (basically, the bary- 
onic physics act as a buffer between the halo-halo merger rate 
and galaxy-galaxy merger rate). But clearly, simply adopting 
the halo trees without accounting for how this processing is 
different from what is typically done (where e.g. halo mass 
ratio is the maximum pre-interaction halo mass ratio, for ex- 
ample) could lead to an apparent deficit of mergers (where the 
same mergers are, in fact, present). 

This is related to another illustration in Figure |4] showing 
the result of adopting what might seem to be the most obvious 
definition of halo-halo mergers given a merger tree. Using a 
pair of closely spaced snapshots, and a descendant halo in the 
later snapshot, we simply take all the halos-i-subhalos in the 
earlier snapshot and define the mass ratio as the ratio of their 
instantaneous mass in that snapshot (the last where they are 
independent) to that of the main (most massive) progenitor. 
Specifically, we apply this algorithm to the standard Mille- 
nium simulation MPA merger trees. M ergers are ag a in dra - 
matically suppressed; this was noted in Bun dy et"an ([2007'). 
The issue is that the definition of "halo-halo merger" is not 
trivial: in the merger tree used to compute this particular ex- 
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ample, halos "survive" in the tree so long as they can still 
be recognized as distinct subhalos/substructure/density peaks. 
In effect, this query is measuring the merger rate in terms of 
the instantaneous mass of a secondary subhalo just before it 
is no longer identifiable as a distinct subhalo/structure in the 
simulation (typically well after the two original halos have 
overlapped and become part of the same FOF group). By this 
time the smaller halo will have been stripped significantly, and 
with infinite resolution this particular way of defining merg- 
ers would make all mergers have a secondary mass — > 0. The 
correct definition of secondary mass is clearly some maxi- 
mum mass or mass before interaction/joining into the same 
FOF group. 

Another, somewhat related caveat, regards the time- 
stepping of simulation outputs. We illustrate one of the possi- 
ble dangers in Figure |5] Obviously, for a process to be well- 
resolved in time, simulation or model outputs must be spaced 
with some Af which is much less than the characteristic evo- 
lutionary timescale of interest. For galaxy halos, this is the 
Hubble time at each z- Halos grow with a characteristic dou- 
bling time that is of order the Hubble time (only weakly de- 
pendent on mass); also, at each redshift, the time between 
mergers is of order the Hubble time; and the timescale from 
dynamical friction for a major merger to complete is a factor 
of 1/10 - 1/5 of the Hubble time. In short, Ar < tniz) 
is required, or in terms of redshift, Az ^ 1 at all z where re- 
solved dynamics are desired. 

Without such resolution, mergers will be entirely missed 
- for example, a system with initial mass nii — m could 
experience a merger with mass mo = 3 m, then the product 
(mass = 4 m) merge with a mass mj, = 12 m system. The 
central system has experienced two consecutive major merg- 
ers (/i = 1/3) in this case. However, if the timestepping is 
large, the outputs of the simulation will only indicate that 
at some initial time, halos of mass mi = m, m2 = 3 m, and 
m^ = 12 m were separate, then at a later time they were all 
together, and the most common assumptions applied in re- 
constructing merger histories will assume each merged onto 
the "primary" (^3) separately, giving two minor mergers (a 
1:12 and 1:4 merger, as opposed to the correct 1:3 and 1:3 
merger). Considering the merger history along all branches 
of a given halo tree, the resulting uncertainty in merger rates 
can grow quickly if the time-stepping used is too large. 

If one adopts a "merger timescale" to follow the galaxies 
in their last stages of infall after the subhalo is unresolved, 
then information is obviously lost without narrow timesteps - 
the merger time will necessarily be based on the last resolved 
dynamics, at large radii, which tends to lead to a significant 
over-estimate of the remaining "merger time" (see the discus- 
sion in § |3.3.3t . This and the above can be significant con- 
cerns at high redshifts in semi-analytic models that interpolate 
between individual snapshots of the dark matter background 
at various times. 

Moreover, if merger rates or mass ratios are defined in terms 
of the "descendant" mass - i.e. if the mass ratio is defined 
by the mass of a progenitor m, (in one output) to the mass 
of the final halo ot/ (in the next output), as is common in 

" These issues are described in detail at 
|http : //www . astro . utoronto . ca/~-bundy/millennium/' 
(K. Bundy). Tliere, the authors both outline the definitions that lead to 
this artificial merger rate suppression and present a modified query of 
the Millenium database (courtesy of V. Springel & S. White) that defines 
halo-halo mergers in terms of the more proper pre-accretion secondary mass, 
that does not suffer from this problem. 



many analyses, or if the halo mass bin for which the merger 
is "counted" (in e.g. asking the merger rate at a given mass) is 
m f - then there can be a significant bias without narrow time- 
stepping because the halo will grow by a non-trivial amount 
in mass between the two timesteps (not just from the merger 
itself). 

Figure |5] illustrates this. We extract the halo-halo merger 
rate from the Millenium simulation using a simple query - 
for every halo in a given timestep, we ask whether it had a 
progenitor (other than its primary /main-branch progenitor) in 
the previous timestep with mass > 1/4 the final mass (corre- 
sponding to a > 1/3 merger). Mass here, to avoid the other 
problems above, is defined as the maximum pre-accretion 
mass. But we vary the length of that timestep. If the timestep 
is large, then at redshifts z> 2, the apparent merger rate ob- 
tained by such a method plummets, even though it is well- 
known that the halo-halo merger rate continues to rise in- 
definitely with redshift. This is because the time-spacing of 
simulation outputs becomes large (Az <C 1 is no longer satis- 
fied) and halos grow non-trivially between timesteps, so ma- 
jor mergers may represent only a small fraction of the mass of 
the "final" halo after a large time (in units of the Hubble time) 
passes. 

3.3. Substructure: Tracking Subhalos and Their Dynamics 
3.3.1. The Subhalo Merger Rate 

Of course, a halo-halo merger is not a galaxy-galaxy 
merger; there will be some finite time where the galaxy is 
part of a satellite/substructure system before the final galaxy- 
galaxy coalescence. To follow this, simulations have made 
great improvements in tracking not just primary halos but also 
substruc ture in those halos, identifying subhalos on the fly 
(see e.g. lSpringel et al.ll2001h . Figure |6]plots the merger rates 
of such subhalos, measured from cosmological dark matter 
simulations: the subhalo "merges" when it is effectively de- 
stroyed in the simulation (it is mixed/stripped to the point 
where no subhalo can be identified by the algorithm, at least 
above the simulation mass resolution limit; we discuss some 
aspects of these limits below). The subhalo is continuously 
losing mass to stripping after the initial halo-halo merger (and 
by definition its mass is effectively zero when it "merges"), 
so the mass ratio is most appropriately defined in terms of the 
ratio of the maximum secondary halo mass (maximum mass 
it had when it was its own halo before any onset of stripping) 
to the instantaneous primary mass (technically primary minus 
subhalo mass, since we define an equal mass merger as 1:1, 
or^haio < !)■ 

We show the results of this determination from several dif- 
ferent sources: the subhalo trees from the Millenium sim- 
ulation (the "MPAhalo" catalog); the comparable but inde- 
pendent dark matter simulations in IStewart et al . (2009a|; 
and fr om a third set of simulations presented in I Wetzel et all 
(|2009|) . We also consider subhalo tracking using a slightly 
different methodology, beginning with e.g. the subhalo mass 
function (rather than the differential subhalo merger rate) in- 
side a radi us < R^„, as de termined from cosmological sim- 
ulations in IKravtsov et al.l 12004) . and evolving these sub- 
halos forward from this already small radius according to 

Implicit in Figure|6]is the assumption that all subhalos merge with the 
central one, not with each other This is not strictly true, but simulations find 
it effects merger rates at only the ~ 10 — 20% level, much less than the dif- 
ferences between different models shown t Wetzel et alji200ft : lAngulo et aH 
12001 . 
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Fig. 6. — As Figure|3] but comparing halo-halo merger rates (from'Fakhouri & Ma"200^ to subhalo-halo merger rates. Subhalo mergers are defined when 
subhalos ai'e fully merged or destroyed with /^halo = ^subhalo /^halo; where M(,aio is the instantaneous total primary halo mass, and M^ubhulo = ^infall is the 
subhalo's infall/maximum pre-accretion (pre-stripping) halo mass. We compare with the results from the subhalo trees in the Millenium simulation (the same 
simulation as Fakhouri & Ma 2008), taken from Angulo et al. (2009); Kitzbichler & White (2008); with the subhalo trees constructed from the independent 
DM simulation in Stewart et al. (2009a); from the subhalo mass functions determined in Kravtsov et al. (2004) with the methodology outlined in Hopki ns et aTl 
t2008b): and with the subhalo mergers defined in another different manner in Wetzel et al. (2009) For the Wetzel et al. (2009) results, shaded range (left) shows 
typical statistical and systematic uncertainties; double lines (right) bracket this range, they are broken as a function of redshift because of the different volume 
simulations used to probe different redshifts. Despite the fact that subhalos can survive a significant time after halo-halo merger, the rates are similar within a 
factor ~ 2 to the halo-halo merger rates; different subhalo identification methodologies also yield similar results to within the same factor 



the results of high-resolution galaxy-galaxy merger simu- 
lations used to calib rate e.g. dynamical friction times (see 
iHopkins et aO2008bllal) . 

Despite the fact that these are all independent simulations, 
with somewhat different algorithms used to identify and track 
subhalos, the agreement is reasonably good: differences in the 
major subhalo merger rate are at the factor < 2 level. More- 
over, the predicted subhalo major merger rate is quite similar 
to the predicted halo-halo major merger rate; depending on 
the calculation, subhalo major mergers appear to be between 
a systematic factor of ^ 2 lower than, or equal to, the halo- 
halo major merger rate. 

3.3.2. The Halo Merger Rate with a "Merger Delay" 

Another common approach, especially given low- 
resolution cosmological simulations or analytic (e.g. 
extended Press-Schechter) approaches where subhalo track- 
ing is not possible, is to assign some "merger timescale." This 
is a delay between the initial halo-halo merger and the final 
galaxy-galaxy merger, approximating the time for the decay 
of the secondary orbit and final resonant tidal interaction 
and coalescence of the two galaxies. Typically some variant 
of the standard Chandrasekhar dynamical friction time is 
employed, but this does not have to be the case: in fact, the 
angular momentum loss and orbital decay in the final stages 
of the merger are governed by resonant tidal processes, not 
dynamical friction against the smooth background. Because 
of these uncertainties, a number of authors have attempted to 
calibrate these timescales using high resolution simulations of 
galaxy-galaxy mergers and halo-subhalo or group encounters: 
we compare such calibrations in Figure]?] 

We consider: (a) The directly extracted time required for 
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Fig. 7. — Comparison of different calibrations of the "merger timescale" 
(average time between halo-halo merger defined at crossing of Syir ™d cen- 
tral galaxy-galaxy merger) in units of the Hubble time as a function of initial 
mass ratio. Each comes from high-resolution N-body simulations: dissipa- 
tionless binary galaxy+halo merger simulations surveying orbital parameters 
from Boylan-Kolchin et al. (2008) (thick line corresponds to the range ex- 
plicitly simulated; thin dotted line to extrapolation beyond this range), hydro- 
dynamic binary galaxy+halo mergers (Cox et al. 2006; Younger et al. 2008}; 
lower-resolution hydrodynamic cosmological simulations ( Jiang et al. 2008!), 
and dissipationless gravitational/group capture cross sections for galaxies in 
groups ( Krivitskv & Kontorovich 1997). All yield a similar factor ^ 2 scatter 
in timescales owing to cosmological distributions of orbital parameters (lines 
are for median circulaiity and other orbital parameters in cosmological sim- 
ulations). Different estimators agree well where calibrated. Major-merger 
timescales are <C fHubble^ as a result, subhalo major merger rates are not very 
different from halo major merger rates. 



the final merger in a galaxy-galaxy hydrodynamic (SPH) 
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merger simulation (of disk+bulge+halo systems), including 
star forma tion, rnulti-component galaxies, and gas physics 
jCox et alJ 120061: 1 Younger et alj 120081: [Hopkins et al.]l2009al; 
lLotzetalJl2008ah . We date the final merger by the coales- 
cence of the two galactic nuclei. Error bars show the variation 
in timescales corresponding to sampling the full (isotropic) 
range of relative initial disk inclinations and orbital parame- 
ters (prograde mergers being a factor ^ 2 more rapid to coa- 
lesce than retrograde mergers, owing to the enhanced tidal re- 
sponses dissipating angular momentum), (b) Calibration of a 
revised dynamical friction formula fro m gas-free (dissipation- 
less) bulge-nhalo merger simulations in iBoylan-Kolchin et alJ 
JToOSi) . These authors also survey a wide set of orbital param- 
eters: we show the result for median and ± 1 cr range of orbital 
parameters measured in cosmological simulations ( BensonI 
E)05; Khochfar & Burkert 2006). We show the full fit, but 
note that the simulations used span /^haio = /^gai = 0.025 — 0.3, 
over which range they agree with the SPH simu lations, (c) 
A sim ilar fit, for cosmolo gical SPH simulations iJiang et aU 
(120081) & lJiang et alJ jloi O). calibrated from /ihaio = 0.04 - 
1 .0. The cosmological nature of the simulations allows full 
sampling of representative orbits, with self-consistent galax- 
ies; however the lower resolution suppresses resonant effects 
in major mergers, contributing to a slightly longer merger 
time, (d) The characteristic timescale for pair-pair gravita- 
tional capture in lo ose group or field environments, c alibrated 
from sirnulatio ns in lKrivitsky & Kontorovichl(ll997h (see also 
iMamonI l200 6l). Although dynamical friction can be a rea- 
sonable approximation to these scalings, such capture cross- 
sections are technically more appropriate for collisions in e.g. 
small groups or field environments, or major mergers, where 
the idea of a long-lived, slow inspiral is not appropriate. An 
alternative cali bration in angular-momentu m space gives sim- 
ilar results (see lBinney & Tremaindll987h . 

These ti mescales are discussed in greater detail in 
iHopkins et al. (2010) andHopkins et al. (2008b). For our pur- 
poses they all agree reasonably well over the range where they 
are calibrated - to a factor better than ~ 2. Where they are 
extrapolated, of course, they should be regarded with some 
caution. 

Figure [S] shows the results of combining these merger 
timescales with halo-halo merger rates. In detail, each halo- 
halo merger is assigned a "merger time" based on the curves 
in Figure [T] appropriate for the redshift of the initial halo- 
halo merger and the halo-halo mass ratio at that moment, as- 
suming the secondary begins at R^h (which is what the fit- 
ted merger times are calibrated to use). The final merger is 
then delayed by this interval, during which the primary halo 
continues to grow and accrete. The mass ratio of the final 
merger is defined as the ratio of the maximum (pre-accretion) 
secondary halo to the primary halo mass at this time of fi- 
nal merger (as is done for subhalo me rger rates). We con- 
sider b oth the fitted merger ti mes from iBovlan-Kolchin et alj 
(l2008h and lJiang etaP (l2008l) . which bracket the other calcu- 
lations in Figure [T] 

To see how merger rates would be affected by very 
long merger timesc a les, w e have also considered the 
iBoylan-Kolchin et alJ (l2008h time multip lied by a constan t 
factor of 4 (for major mergers, roug hly, the lJiangetal] (l2008h 
time multiplied by a factor of ~ 2, or '-^ 0.5fHubbie(z))- This 
represents a fairly extreme case. We compare to the subhalo 
and halo-halo merger rates, as a function of mass ratio and 
redshift. 

The merger rates calculated with this approach agree to 



within a factor ^ 2 with those determined from a direct track- 
ing of subhalos in simulations: it appears that evolution on the 
timescales adopted for "typical" orbital parameters in simula- 
tions, shown in Figure |7] is a good approximation to the full 
dynamics. The redshift evolution is slightly different: high 
redshift mergers (when halo growth is more rapid) have their 
mass ratios somewhat more suppressed, but in terms of merg- 
ers per unit redshift the evolution is still weak and within the 
range of various subhalo calculations.'^ The differences be- 
tween each different merger time calibration are small (for 
major mergers): it is only when we artificially increase the 
merger timescale by a factor of several that a significant dif- 
ference appears, in the sense that high-redshift mergers are 
strongly suppressed. This emphasizes that, although exist- 
ing calibrations from high-resolution simulations agree well 
and yield little difference when used properly, adopting a sig- 
nificantly longer merger time can substantially suppress the 
merger rate, especially at high redshifts. 

3.3.3. Caveats: Subhalo Identification and The Application 
of Dynamical Friction to Orphans 

Despite the reassuring agreement seen above, it should be 
noted that merger rates are sensitive at a non-trivial level to 
issues of subhalo identification: in particular the choice of 
algorithm and criteria used to determine subhalo "member- 
ship," and the time, mass, and force resolution of the sim- 
ulation. These effe cts have been examined in a num ber of 
simulations (see e.g. [K lvpin et al. 1999; Springel et al. 2001; 
IShaw et alllioOTl: IWetzel et al. 2009; W etzel & White. 20101: 
Gioco li et al.l]20ia " and references therein). We note here a 
couple of possible pathological regimes. 

First, it is possible in some particular circumstances that 
subhalos could survive "too long" in the simulations shown 
above, depending on the definition of the time of merger. 
This can happen, for example, because the contribution of 
baryonic galaxies to speeding up the merger is neglected. 
Including baryons keeps subhalos more tightly bound and 
changes the mass distribution, shortening the dynamical fric- 
tion time (see Weinberg et al. 2008). Moreover, the baryonic 
components actually control the final merger, allowing for 
resonant interactions once the secondary is within a radius 
that contains a mass less than or of order its own mass (at 
which point the Chandrasekhar dynamical friction approxi- 
mation is no longer valid). This makes the merger time within 
^ 20 — 50kpc shorter by a factor « 2, re lative to that esti- 
mated without baryonic galaxies (see e.g. iLotz et al.ll20b8al 
2010b a), and can even lead to the merger and/or destruction 
of the baryonic galaxies before the subhalo (their bound dark 
matter) i s significantly g ravitationally disrupted or stripped 
(see e.g. iD'Onghia et all 12009). The lack of very high res- 
olution and accurate ab initio formation of realistic disks (in 
e.g. their scale heights and sizes) in full cosmological simula- 
tions contributes to the fact that the major merger timescales 
in the calibration of Jiang et al. ( 2008) are significantly longer 
than those in high-resolution hydrodynamic simulations of 
individual mergers of disk galaxies. And in the opposite, 
infinite-resolution limit of dark-matter simulations, a subhalo 

" If satellite-central mergers were simply a delayed version of halo-halo 
mergers, keeping mass ratio fixed in infall, then the merger rate distribution 
versus mass ratio would have the same shape as that of halo-halo mergers but 
with a redshift lag - and because merger rates increase with redshift, it would 
be higher than the immediate halo-halo merger rate. But differential mass 
evolution in this time leads to fewer major mer gers. Together these drive the 
effects seen in Figure [SlCsee lWetzel et alJ2009t) . 



14 Hopkins et al. 




M-haio [foi" SubHalos] z 

Fi g. 8. — C omparison of halo and subhalo merger rates from Figure |6] to the merger rate determined by taking an "initial" halo-halo merger rate (that 
from IFakhour i & Ma 2008) and combining with the merger timescales from Figure [7] (i-e. including a "delay" between the time of halo-halo merger in the 
simulation and the final recording of the merger or assumed subhalo-halo merger, given by the different merger timescale fits in Figure [T)- The timescales of 
[Boylan-Kolchin et al. (2008) and Jiang et al. (2008) bracket the major-merger timescales therein. The difference is relatively small. We also compare the result 
given by artificially increasing the merger time by a large factor (= 4) relative to the calibration o f Bovlan-Kolchin et aL ti2008i) : in this case the merger timescale 
becomes so long (r^ fHubble) that high-redshift mergers are substantially suppressed. 



might never be completely disrupted (there always being a 
few tightly bound particles), even when the contained galax- 
ies would have long ago merged. 

Correcting for these effects is difficult. Subhalos could 
be identified as "merged" if they hav e lost some large frac- 
tion (> 90 - 99%; for calibrations see I Wetzel & Whitell2010l) 
of their mass (as opposed to simply reaching some resolu- 
tion limit), but in low-resolution simulations this may prema- 
turely "destroy" the subhalo. One could attempt to account 
for the resonant effects of baryons by merging any systems 
that have a passage within a close radius r to another galaxy, 
such that the enclosed mass M{< r) is comparable to the sub- 
halo/secondary mass (really the secondary galaxy mass plus 
tightly bound dark matter within ^a couple R^) - within any 
such radius the resonant and baryonic interactions will be effi- 
cient and dominate over the dark matter dynamics and simple 
dynamical friction is not applicable. 

A more common problem, especially in cosmological sim- 
ulations (which by necessity have less than ideal resolution), 
is that subhalos are lost below the resolution limit "too early." 
Identifying a subhalo as a distinct bound substructure is de- 
manding in terms of resolution, and at some point subhalos 
will fall below these limits as they are stripped. Lack of inclu- 
sion of baryons (which would keep the subhalo more tightly 
bound), lower force resolution, or conservative subhalo mem- 
bership criteria will all accelerate this. Subhalos in the Mil- 
lenium simulation, for example, typically fall below the res- 
olut ion limit within ~ 0.5 — l^vir (at low masses in particu- 
lar). I Wang et aP (l2006h show that if one simply took only the 
still-resolved subhalos in this simulation at any given time, it 
would not be possible to reproduce the clustering/halo occu- 
pation of observed galaxies on small scales (< 100 kpc), or 
that of subhalos in higher for ce-resolution simulation s. Simi- 
lar conclusions are r cached in I Wetzel & White! (|2010|) . 



To compensate for this, many models built on halo-i-subhalo 
trees invoke a population of "orphan" galaxies. This is the 
same principle as the "merger timescale" above: the models 
follow subhalos in the cosmological simulation as far as possi- 
ble, and then assign a "merger time" to the galaxy based on the 
last recorded subhalo properties in the last timestep where the 
subhalo could be identified. Usually, the assigned timescale 
follows the Chandrasekhar dynamical friction formula. 



^merger 
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I'vir'sat 



Gm.satln(l +Mvir/msat) 



= (1.05(e)"*' +0.6) 



(1) 

1/2 



Vvir ln(l +Mvii-/msat) \^v 

from iBinnev & Tremain3 ( Il987h, with the second equality 
coming from the calibration in I Jiang et aTl (l2008h assuming 
initial rsat = ^vir (their Equation 8; e = j/ jc is the orbital ec- 
centricity with (e) = 0.5, and Rc/R^u is related to the initial 
pericentric passage distance, both taken from the distributions 
they measure). Fits to high-resolution simulations with some- 
what more freedom in the scalings give 



(Mvir/msat) 



1.3 



(—^ (2) 
i(l+M™/msat) V/?v„-/ 

from ' Bovlan-Kolchin et alj (120081) (again adopting a the me- 
dian initial circularity j/ jc = 0.5), where m^sx and rsat are the 
satellite mass and location at the last time its subhalo could 
be identified, and Vvii and Mvir are the virial velocity and mass 
of the primary (minus the subhalo/secondary itself). Similar 
scalings are obtained from analytic considerations (at least in 
the small msat/A^vir limit), accounti ng for subhalo mass l oss 
and orbital parameter evolution (see Taylor & Babul 2004). 

These particular formulations, discussed above (see Fig- 
ure |7]), have been tested in simulations and shown to be a 
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Fig. 9. — Caveats of using a subhalo+merger time ("orphan") calculation of 
the total merger timescale. Top: Ratio of the "subhalo plus oiphan" timescale 
?nigi (sub + orphan) to the true merger timescale f,„gi (full simulation) from 
Svii to final merger The true timescale fmgr (full simulation) is given by the 
calibrations in Figure [T] - simulations that follow galaxies fully from Ry„ 
to R — > 0. The subhalo+orphan timescale comes from following a subhalo 
self-consistently to some i?(subhalo destmction) where the subhalo is lost 
(falls below some resolution limit) - a timescale given by the same cali- 
brations - then re-applying the merger time formula based on the instan- 
taneous S(subhalo destruction) and remaining galaxy-l-subhalo mass. The 
ratio of this timescale to the true merger time is shown as a function of 
R(subhalo destruction), the initial halo-halo mass ratio, and galaxy mass. 
If R(subhalo destruction)/Svir > 0.2, the total merger time can be signifi- 
cantly over-estimated. Middle: Average ^(subhalo destruction) relative to 
i?vii , where subhalos in the Millenium simulation fall below the mass resolu- 
tion limit, as a function of the secondary/satellite halo infall mass and redshift 
(S. Genel; private communication). Small satellites, which are close to the 
simulation resolution limit, naturally fall below the resolution hmits at rel- 
atively large R in any simulation. Bottom: Effect on the "typical" inferred 
merger time for 1 : 3 mergers, assuming subhalos are lost at the radii of the 
z = curve at middle and given the methodology shown at top, as a function 
of galaxy mass (using the Mgj,i — Mi,aio relation from Wang et al. 2006). If 
subhalos are only marginally resolved, caution should be used in applying 
merger times based on instantaneous subhalo properties; better calibration is 
needed. 



reasonable approximation to the total merger time from ini- 
tial halo-halo contact at rsat — 1 — 2/?vir, with Wsat defined 
as the maximum (pre-stripping) mass of the secondary (this 
is the input when used to calculate the merger time for a 
halo-halo merger, as in § |3.3.2| i. However, these formula- 
tions have not been calibrated for application to a satellite at 
a later time, after the subhalo has self-consistently evolved 



and is heavily stripped. Assume, for example, that a subhalo 
with initial (pre-stripping) mass m, self-consistently decays in 
orbit from an initial radius Rj ^ Rvk acc ording to the calibra- 
tion from B ovlan-Kolchin et~an ( l2008l) (Equation |2] above), 
but at some point - when the subhalo is at a radius r'^ = 
7? [subhalo destruction] and time t — t^^^ |2](^' ^ '"s I ~ '"0 
- can no longer be resolved. At this point the subhalo mass 
is mj. Now assume that the remaining merger time is then as- 
signed according to the same Equation |2l with these revised 
input values, i.e. a remaining t — f^^^ ^''s -^0\m — m'^). 

The total merger time is then simply the time from /?vir to 
plus the assigned remaining time. But we also know what the 
total merger time should be from the initial /?vii radius to r = 
(i.e. without breaking into sub-steps in this manner) - this is 
in fact what the formulae used above are explicitly calibrated 
from high-resolution simulations to reproduce (with the in- 
puts of initial /?vir and m,). Comparing the two, we obtain the 
ratio of the total merger time estimated from this two-stage 
breakdown, to the total merger time from 7?vir to r = directly 
from the simulations. This is 
?mgr(sub + orphan) 



?mgr(full simulation) 



fEqn.|2](^vir "> | '« = '«,) + t^^^^ ^r', -^0\f71 = m',) 



t^ 



1 



'Eqn.|2](^vir^0|m 

il^ ln[l+Mvir/m,] 



;<i-3ln[l+MvirM] 



But nis/Myir is just the initial halo-halo mass ratio (/i,). The 
final mass m'^, evaluated instantaneously when the subhalo is 
nearly or entirely stripped, will by definition just be the galaxy 
mass Mgai = ?y//,Mhaio, where fh is the Universal baryon frac- 
tion, Mhaio is the pre-accretion subhalo mass, and ?/ is some 
star formation efficiency (77 = 0.3 around ^ where star for- 
mation is most efficient and drops rapidly at lower or higher 
masses). If we include the mass of the subhalo at the last time 
it was resolved, we add to this roughly the simulation resolu- 
tion limit; but for intermediate to high mass galaxies and and 
the resolution limits of simulations of interest, this is similar 
(^ 10"^ — IO'^Mq). For typical ^ L* galaxies and the parame- 
ters of e.g. the Millenium simulation, with an initial 1:3 halo- 
halo merger, this yields f(sub + orphan)/? (full simulation) « 
1 + (6.5 X — l)x, where x — r'jR^^y is the fraction of the virial 
radius where the subhalo is no longer resolved. 

Figure |9] illustrates this as a function of initial mass ratio, 
galaxy mass/resolution threshold, and the radius where the 
subhalo becomes unresolved. In short, if subhalos can be self- 
consistently followed to small radii < 0.1 — 0.2/?vii, where the 
remaining merger time is short, there is no issue. However, 
if subhalos fall below the resolution limit at large radii, then 
the merger time can be over-estimated by a large factor. We 
compare the typical radii in the Millenium simulation where 
subhalos of a given initial (pre-accretion) mass become un- 
resolved (after being accreted into a larger halo), at several 
redshifts (S. Genel, private communication). Unsurprisingly, 
when the subhalo initial mass is close to the resolution limit, 
subhalos are lost at ^ R^h. Given the Mgai(Mhaio) relation and 
these median radii as a function of mass, the Figure illustrates 
how this affects the expected merger time for satellites of a 
given mass: at high masses, subhalos are massive and can be 
followed self-consistently, at low masses however, subhalos 
are lost quickly and the application of Equations [T]|2] to the 
"orphan" tends to over-estimate merger timescales. This will 
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have a large effect on predicted galaxy mass assembly at high 
redshifts. 

The issue is that fitting formulae such as Equations [T] & 
|2] are not necessarily calibrated for these multi-tiered appli- 
cations. This is not to say that dynamical friction is fun- 
damentally invalid - orbits of small satellites can be well- 
approximated by application of an insta ntaneous dynamical 
friction drag force and tidal stripping touinn & Goodman' 
1986; Hernquist & Weinberg 1989; Velazquez & White 1999; 
Font etal.1 1200 It jVillalobos & Hehnii l2008t iBenson et alj 
2004|). (Although for major mergers, the underpinning as- 
sumptions of dynamical friction are not strictly valid). But 
this is not as simple as a single timescale, and implicitly as- 
sumes a high-resolution limit - they are not necessarily de- 
signed to account for a resolution "floor" self-consistently. 
Moreover, many estimates of the merger time average out or 
ignore the detailed orbital information; often, cases where the 
subhalo is quickly dissociated/stripped correspond to highly 
radial orbits, and simply re-applying an orbit-averaged for- 
mula is not ideal. If a system has just had a radial passage, for 
example, it will most likely be "caught" in a snapshot at apoc- 
enter - applying the dynamical friction time from this point, 
assuming the orbit is circular (as in most models) will in fact 
over-estimate the remaining merger time by a factor of « 6. 

4. THE LIGHT SIDE: HOW GALAXIES OCCUPY HALOS 

The issues discussed thus far pertain to all models, whether 
they are halo occupation, semi-analytic, or simulation-based. 
Although we have seen that some of the aspects of the treat- 
ment of dark matter mergers can be important, we know from 
our comparison in § |2] and the discussion above that most 
of the order-of-magnitude differences in merger rates must 
primarily owe to the treatment of baryons. Once we know 
how halos/subhalos behave, then the other necessary ingre- 
dient to predict galaxy merger rates is some prescription for 
how galaxies populate those halos (in other words, the "halo 
occupation" statistics). We therefore investigate how differ- 
ences in these halo occupation aspects of the models give rise 
to different merger rates. 

4.1. Halo Occupation Statistics of Central Galaxies 
4.1.1. At Low-Redshift 

With some complete a priori knowledge of subhalo evolu- 
tion and mergers, the galaxy-galaxy merger rate, as a function 
of galaxy mass and mass ratio, can be determined knowing the 
distribution of galaxy masses in each halo. To lowest order, 
this is given by the function Mgai (Mhaio , z) and its scatter. Of 
course, this could depend on other variables such as environ- 
ment, but as discussed below, observations indicate that such 
additional parameters have much smaller systematic effects 
than the relevant uncertainties in Mgai (Mhaio , z) ■ 

For simplicity, and to isolate the important physics, we 
first consider only central galaxies (i.e. those where Mhaio 
is the parent halo, not satellites in subhalos), at z = 0. Fig- 
ure[TO] shows the median Mgai(Mhaio), as a function of Mhaio, 
at z = (we could also compare the scatter in this quan- 
tity, but all the models we consider predict a similar, rela- 
tively small 5, 0- 15 — 0.2 dex scatter; see e.g .lYang et al. 2008; 
Mosteretal.' '2010b'; 'Behroozi et al.' '2010; TWetzel & White! 
2010; More et al. 2009; Guo et al. 2010^^ Of course, since 
every satellite was a central galaxy in some halo of mass 
Mhaio before being accreted, the halo occupation function and 
Mgai (Mhaio) function of satellites should be identical (barring 



strong satellite-specific physics and/or strong redshift evolu- 
tion in the Mgai (Mhaio) function since the time the satellite was 
accreted), so long as Mhaio for satellites is defined as the max- 
imum (pre-accretion) mass of the satellite subhalo (i.e. max- 
imum mass when the satellite was central in its own halo). 
We w ill free this assumption and discuss the implications in 
§ 14.21 below. It is worth noting that the vast majority of the 
galaxy number density and mass density at all galaxy masses 
is in central galaxies: this is the quantity that is in practice 
most strongly constrained by a model's fit to the galaxy stel- 
lar mass function or the clustering amplitude (large-scale bias) 
of galaxies as a function of mass/luminosity. 

Given some function Mgai (Mhaio), the galaxy-galaxy merger 
rate is simply this function convolved with the population 
of halo-halo mergers, with the appropriate accounting for 
delays before the final subhalo merger/destruction or the 
"merger timescale" from e.g. dynamical friction. Since we 
aim to isolate the specific effects of differences in the distri- 
bution Mgai (Mhaio), we show in Figure[TO]the resulting merger 
rate after convolving each presented halo occupation func- 
tion with th e same halo merger rate (in this case that from 
lFakhouri~& Ma 2008 ) . A different choice of halo merger rates 
or subhalo treatment will systematically change the predicted 
rates along the lines shown in §[3] but will not change how the 
galaxy-galaxy merger rate systematically depends on HOD 
quantities. 

We compare the predicted Mgai (Mhaio) distributions from 
different theoretical models with several observational esti- 
mates. First, fitting to the correlation function of galaxies with 
a given baryonic or stellar mass^" yields a characteristic host 
halo mass for each galaxy mass bin - we show the results 
of this exercise p erformed on SDSS galaxies at z = from 
IWang et aP (l2006l) . We also show the empirical "abundance 
matching" or "rank ordering" results: it has been shown that 
good fits to halo occupation statistics (e.g. group counts, cor- 
relation functions as a function of galaxy mass and redshift, 
etc.) over a range of redshifts z — 4 are obtained by sim- 
ply rank-ordering all galaxies and halos-i-subhalos in a given 
volume and assigning o ne to another in a monotonic one- 
to-on e manner (see e.g. Conroy et al.l 12006'; 'Vale & Ostrike^ 
l2006l) . Here we perform this exercise using the redshift- 
dependent stellar mass functions from Fontanaet al. (200^ 
(showing the result at z = 0); but as discussed in § 12.21 
other choices of stell ar mass function make little differ- 
ence. Recent studies (Guo et al."2009^. 'Mos ter et alj|2010bl; 
[Tinker & Wetzel 2010; Behroozi et al. 2010) obtain similar 
results (they are statistically indistinguishable from the other 
observational results shown in Figure [Toli. We compare these 
results to independent estimates of the Mgai (Mhaio) distribu- 
tion, whe re halo masses a re determined from weak lensing 
studies in lMandelbaum et al. (2006). Other independent mea- 
surements give nearly identical constraints: these include es- 
timates of halo masses in low-mass spirals as a function of 
mass from rotation curve fitting (see e.g. iPersic & Saluccil 

Because it is closer to the relevant parameter in e.g. merger simula- 
tions, and because it better highlights the relevant systematic effects, we 
will show our comparisons here usually in terms of the galaxy baryonic 
massMg;,]. However, our qualitative statements in this section apply to stel- 
lar mass M, as well. Where necessary, we will convert between observed 
stellar and baryonic masses using the fits to the observed Mgas / M, relation 
of g alaxies as a function of redshift from z = — 3 given in Stewart et al] 
12009b); this is consistent with a large number of direct observations over 
this redshift range (Bell & de Jons 2000; McGauah 2005; Calura et al. 200|; 
Shapley et al, 2005 ; Erb et al. 2006; Puech et al . 2010; Miannucci et al.2009l: 
ICresci et alJ2009t [lister Schreiber''enill)09l;lErbl2068l) .'' 
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log(M,„„/h-^M<,) logMg„[Md 

Fig . 10. — Left: Halo occupation statistics (mean galaxy baiyonic mass as a function of host halo mass) from observations of weak lensing ('Mandelbaum et ajj 
120061) . fits to halo occupation models from clustering data in Wang et al. (2006), using the monotonic abundance matching method in Conroy et al. (2006), or 
from group catalogs with satellite abundances and kinematics (Yan g et aL.2008) . We compare with different semi-analytic model predictions, and the simplest 
efficient star formation model (Mgjii = fbM\^.^\„; similar to hydrodynamic simulations without feedback). Right: Corresponding merger rates (at z = 0) as a 
function of M„ai (in ±0.5 dex bins of A/j^i). combining these HODs with the dark matter merger rates from Fakhouri & Ma ( 2008). We hold all other properties 
fixed. For now, we consider just the central galaxy HOD, and assume satellites obey the same relation (where Mhaio is the maximum pre-accretion mass for 
satellites). Different observational constraints yield small (factor ~ 1.5 near ~ L») differences at z = 0; SAMs (being constrained to reproduce the stellar mass 
function) agree within factors ~ 1.5 — 2. However, simply taking halo-halo merger rates or A/j^i = /fcA/haio, the merger rate predicted is different by factors of 
~3-10. 



19881 [19901: iPersic et al.ll996bllaL lBell & de Jond200dl200Tt 
Borriello & Salucci 2001: Borriello et al. 2003l: IShankar etall 
2006: Avila-Reese et al.. ,2008 ) and estimates of the halo 
masses of massive groups and clusters from e.g. X-ray 
gas and group kinematics, as a function of central galaxy 
mass in massive brightest g roup or brightes t cluster galax- 
ies (BGGs/BCGs; see lEke et all 12 004: Yang et al.1 l2005at 
iBrough et al .1120061 Ivan den Bosch et al.ii2007l and references 
therein). 

In semi-empirical or halo occupation-based models of 
galaxy-galaxy merger rates, the required Mgai(Mhaio) distri- 
bution is simply adopted from these or similar empirical con- 
straints, and enforced on all the galaxies in the model. Such 
models are therefore different only insofar as different em- 
pirical constraints yield different estimates of Mgai(Mhaio)- 
We find that the effects of these differences are relatively 
small (factor < 2 near ^ L*, at z < 2); the level of vari- 
ation within such models ow ing to observa t ional uncertain- 
ties is discus s ed in d etail in iHopkins et al.l (1201 Ol) and also 
IStewart et al.l (l2009ah . and summarized in § 12.21 Figure [TO] 
also clearly illustrates that the different observational esti- 
mates of the Mgai(Mhaio) distribution yield similar merger 
rates; the uncertainties are largest at very low mass (where 
observations are difficult and knowing the true gas fractions 
of galaxies, which are gas-dominated at these masses, is crit- 
ical) and very high mass (where statistics are few and, given 
the stellar and halo mass function shapes, the merger rate will 
be most sensitive to small changes in the HOD). 

Now compare these results with those from the Mgai(Mhaio) 
distribution predicted by various a priori models. The semi- 
analytic models, for the most part, match the observational 
constraints quite well. This should not be entirely surpris- 



ing: a major constraint used in building those models is that 
they reproduce the observed z = galaxy luminosity and stel- 
lar mass functions. Such a result (an important achievement 
in and of itself, we stress) implicitly guarantees a reason- 
able match to the observations in Figure [TO] Of the specific 
choices shown, the Bower et al. (2006) model appears to de- 
viate from the observations to a slightly greater extent, around 
^ L*. The manner of this discrepancy reflects a common ten- 
dency for models to not reproduce quite as sharp a break in 
the stellar mass function as is observed. In a very narrow stel- 
lar mass bin, this could lead to non-trivial differences in the 
predicted merger rate, but averaged over a reasonable base- 
line (~ 0.5 — 1 dex in stellar mass), the shape of the predicted 
HOD is nevertheless in sufficiently good agreement with the 
observational constraints that it yields merger rates converged 
at the ^ 50% level. It seems clear that, insofar as they succeed 
in reproducing the observed stellar mass functions, the differ- 
ences between merger rates in semi-analytic models owing to 
differences in the predicted Mgai(Mhaio) distribution of cen- 
tral galaxies at low redshifts are no larger than the differences 
that result when one attempts to directly adopt the Mgai(Mhaio) 
distribution from observatio nal constraints. 

As discussed at length in iHopkins et"an (1201 Oh . the HOD 
shape is non-trivial, and this leads to a significant galaxy-mass 
dependence of the galaxy-galaxy major merger rate, even 
though the halo-halo merger rate depends only weakly on halo 
mass. If the galaxy formation efficiency were to be indepen- 
dent of halo mass, so that Mgai oc Mhaio, then galaxy-galaxy 
mergers would be a direct, trivial reflection of halo-halo merg- 
ers (regardless of the absolute value of that constant of propor- 
tionality). However, at low masses, the dependence of Mgai 
onMhaio is steep, Mgai ^^1^^^^'^. So a "typical" 1:3 halo-halo 
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merger is a 1:9 galaxy-galaxy merger, and the galaxy-galaxy 
major merger rate is correspondingly suppressed (relative to 
the halo-halo merger rate). At high masses, the dependence 
is shallow, Mgai oc mJJ^j^"''', so a typical 1:9 halo-halo merger 
is a 1:3 galaxy-galaxy merger, and the galaxy-galaxy merger 
rate is enhanced. 

This is important for model merger rates in hydrodynamic 
galaxy formation simulations: it is well-established that in 
cosmological hydrodynamic simulations (at least those with- 
out detailed prescriptions for both stellar and AGN feed- 
back), ga laxy formation is efficient, with Mgai ^ /bMhaio 
(see e.g. ISpringel & Hernqulsg l2003bl: iMaller et all |20^ 
iNaab et aTl I2007L and references therein). Although con- 
siderable progress is being made (see e.g . Governato et all 
2007, 2009; Siiacki et al. 2007; Scannapieco e t all 120(31 
Pi Matteo et al. 2008; Croft et al. 2009 ). no cosmological hy- 
drodynamic simulation has yet been successful at reproduc- 
ing the z = stellar mass function and HOD; and simulations 
without feedback tend to produce something much closer to 
efficient cooling/galaxy formation at all masses, with Mgai oc 
Mhnio (or at least closer to this limit than to the Mgai(Mhaio) 
distribution observed). As a consequence, such a simulation 
will predict a galaxy merger history that simply reflects the 
halo growth history: fewer major (and more minor) mergers, 
with a factor of several bias in the predicted merger rate. This 
effect is especially important at masses > L, ; a cosmological 
hydrodynamic simulation that does not correctly reproduce 
galaxy "quenching" and match the observed bright end of the 
galaxy stellar mass function can under-predict the number of 
major mergers by factors of ~ 5 — 10. 

In fact, these effects have been se en in such simulations 
(IMaller et al.ll2006t iNaab et aLllIOOTl) . Our comparison here 
demonstrates that indeed this is expected in a simulation of 
this nature and reflects halo growth. However, if we take these 
same simulations, and simply re-populate all of the galaxies 
according to the observed HOD constraints (i.e. take the simu- 
lation but replace the galaxies just before merger with ones of 
the "correct" mass according to the observational constraints), 
then we obtain the results from Figure [T] These results agree 
to within a factor ~ 2 with the semi-empirical models and ob- 
served merger rates. 

4.1.2. High-Redshift Results 

Thus far we have only considered results at z = 0. FigurefTTI 
compares the Mgai(Mhaio) distribution as a function of red- 
shift, from different observational estimates and semi-analytic 
models. The agreement is reasonable around ^ L*, but less 
so at higher and lower masses. This reflects well-known as- 
pects of the comparison between such models and observed 
galaxy stellar mass functions - the SAMs tend to predict 
more low-mass and fewer high-mass galaxies at high red- 
shifts (Fontanaet al. 2006; ll bert et al.ll2010t iFontanot et"an 
120091: [Marchesini et al. 2009) . Of course, the uncertainties in 
SAMs grow at high redshift, but so do observational uncer- 
tainties. Low-mass galaxies are subject to completeness con- 
cerns and selection effects; high-mass galaxies may have non- 
trivial corrections from extended intra-group light and/or un- 
certainties in their pho tometric ally inferr ed stellar masses an d 
photometric redshifts (iMarasto n 2005; Maraston et al.ll2006l) . 
and there may be evolution in the stellar initial mass function 
(glopkins & Beacom 2006; van Dokkum 2008; Dave 2008). 

The growing differences and uncertainties in Mgai(Mhaio) 
highlight the need for further study and observational con- 
straints. Nevertheless, the resulting merger rates, integrated 



over a reasonably broad stellar mass interval, are still similar 
to within a factor ^ 2. This owes to the fact that, at a given 
stellar mass M», the merger rates do not depend on the abso- 
lute value of Mgai(Mhaio), but rather the shape of this function. 

4.2. Halo Occupation Statistics of Satellite Galaxies 

We now generalize our previous comparison to allow for 
different Mgn\{Mha\o) distributions in satellite or central galax- 
ies. In other words, whereas before, a satellite galaxy's prop- 
erties were set by its maximum pre-accretion halo mass (i.e. 
its state when it was last a central galaxy), we now free that 
assumption. 

4.2. 1. The Central-Satellite Difference and its Effect on 
Merger Rates 

Figure [12] shows the relevant comparison. Specifically, we 
compare the median Mgai(Mhaio) relation, for both central and 
satellite galaxies. As always, we define Mhaio for satellites 
as the maximum Mhaio before their accretion, i.e. their "in- 
fall" mass Minfau. We sho w two represe ntative models: first, 
the HOD constraints from iConrov et al.l C2006). inferred from 
the observational constraints, in which the two Mgai (Mhaio) 
relations (that for satellites and centrals) are nea rly identi- 
cal (with this pre-accretion defi nition of Mhaio; see iGao et al.l 
2004; Naga i & Kravtsovl 12005). Indeed, that the two distri- 
butions can be nearly the same is one of the major findings 
of such empirical HOD models. To be specific, the obser- 
vational input in this method typically constrains the stellar 
mass M,.sat/M,.cen, not the baryonic mass Mgai,sat/Mgai,cen - 
we assume here similar /gas(M,) for the cold, bound, disk gas 
component. However we will show below (Figure [T3] l that 
observations show no difference in /gas(M*) for the specific 
(major-merger candidate) systems of interest here. Further 
exploration of these distinctions is important, however (e.g. 
Neistein et al., in prep) - for now, this should more conser- 
vatively be considered an assumption in the semi-empirical 
models. 

Observationally, the large-scale clustering of galaxies (and, 
to first order, the stellar mass function) constrain primar- 
ily central galaxies. However, the small-scale galaxy cor- 
relation function (the "one halo" term) probes satellite pop- 
ulations, and independent constraints are available from 
e.g. group counting statistics and satellite kinematics. As 
noted in § 14.11 most halo occupation models find good 
agreement with these observations by assuming that satel- 
lites/subhalos are described by the same Mhaio distribution 
as central galaxies, so long as the subhalo mass Mhaio used 
to assign the subhalo galaxies (satellites) their masses is 
the maximum pre-accretion mass, as oppose d to their in- 
stantaneous (stripped) subhalo mass (see e.g. Conrov et ajj 
2Q06; Zheng et al. 2007; Conrov et al. 2007; Vale & Ostrike] 
2006tlConro v & Wechsler 2009; Perez-G onzalez et al.l2008"al: 
Stewart et al., ,2009a; Hopkins et al. ,201(il). This includes all 
satellites in the model, even those about to merge (at very 
small scales < 20 — lOOkpc), and results in excellent agree- 
ment with the observed number of pair s on those scales as 
a function of ga l axy mass and redsh ift dConrov et al.l l2006t 
iBeU et al.ll200^IStewart et al.ll2009ah . 

The top panel of Figure [12] shows this result: Mgai (Mhaio) 
is identical for both central and satellite galaxies, defined 
in these terms. But although this is a typical assump- 
tion in halo occupation models, the observational constraints 
may not necessarily require identical satellite and central 
properties at fixed infaU mass. Therefore, in terms of the 
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Fig . 11. — As Figure riQl but comparing empirical HOD constraints and a priori model predictions at different redshifts. Top Left: Stellar HOD (median central 
stellar mass versus Mhaio) at z = 1 (lef t) and z = 2 (right). We compare the empirical HOD from abundance-matching, using stellar mass functions/HOD fits from 
[Ojnrov & Wechsler 1 2009); Fontana e t alj 120061) and Perez-Gonzalez et al. 1 2008b). We contrast predictions from the Bower et al. 1 2006), de Lucia & BlaizoJ 
(2007), and Khochf ar & Silki (2009) SAMs. Bottom Left: Same, in terms of baryonic mass. Right: Resulting galaxy-galaxy major merger rate as a function 
of redshift, in a given baryonic mass bin, from the different HOD curves. In general, where the shape of the HOD is similar (regardless of normalization), the 
merger rates are in agreement. 



ratio Msateiiite/ZWcentrai - the typical satellite versus central 
galaxy mass at fixed halo infall mass, the top panel shows 
the allowed scatter or systematic offset that could be toler- 
ated within observational errors , given the constraints pre- 
sent ed in IConroy et al.l (l2006h . IConrov & Wechsled (l2009h 
and IStewart et al. (l2009alK The authors find that less than 
~ 0.2 dex in scatter or systematic offset of satellite popula- 
tions from central populations is acceptable, before agreement 
with the observed clustering o n small scales breaks down (see 
alsolMosteretaLllIO lOb: Gu o et al.1 l2009UTinker & Wetzell 
l2010tlPasquali et al.ir2010i:iWeinmann et al.ll2009bh . 



The conclusion is not strongly redshift-dependent at least 
for redshifts z = — 4 where clustering and stellar mass func- 
tion measurements are available (and observations of other 
populations even at z^6 suggest it remains true the re as well; 
IWhite et al.ll2"008t IShankaJl2009t IShen et al."2010). This has 
also been found in a number of oth er studies (Zhena et al. 
2005, 2007; Yang et al."2008l; ll^ret al.ll2010t IMore et all 
2009; Mo ster et al. .2010bi). 



For example, IWang et aH (l2006l) arrive at similar conclu- 
sions with independent data sets, simulations, and different 
fitting methodologies. They find reasonable agreement (espe- 
cially at low masses, which is the regime of particular interest 
here) assuming that centrals and satellites obey an identical 
Mhaio distribution. They do obtain a marginal improvement 
for high-mass galaxies if the preferred satellite Mgai(Mhaio) is 
offset from the central Mgai(Mhaio) by a small amount, shown 
in Figuredll satellites might be systematically less massive at 
the same halo infall mass by a factor ^1.2 around ^ L^, and 
^2 at L^,. But differences in the fitted scatter actually 



"cancel out" much of the predicted offset, and the statistical 
difference at low masses is small. ^' 

Interestingly, cosmological hydrodynamic simulations pre- 
dict a nearly identical Mgai(Mhaio) for satellite and central 
galaxies, in agreement with the HOD model fits (including the 
assumption that the /gas (M* ) relation is unchanged for satel- 
lites and centrals, and therefore that there is no difference ei- 
ther in M*(Mhaio) orMgai(/Whaio))- The reasons for this are 
discussed below (§ l4.2.3l l. 

We compare these empirical HOD cons traints with the pre- 
dictions of the semi-analytic model from lde Lucia & BlaizotI 
(2007) (the results are similar for other SAMs shown, we 
show this as a representative case). Since we are interested 
in effects on mergers, we specifically isolate populations in 
the "orphan" stage, i.e. for which the "merger clock" in the 
model has been initialized and for which the remaining time 
until merger is less than the Hubble time.^^ In this model, 
and in many other SAMs, there is actually a significant dif- 
ference between the two Mgai(Mhaio) relations (especially at 



^' Note that the simulations used in IWang et al.l (20061) require the inclu- 
sion of an "orphan" (unresolved subhalo but un-merged) population in or- 
der to match the observed small-scale clustering: repeating the analysis with 
subhalo populations in some other simulations (e.g. Stewart et al. 2009a) that 
require a smaller orphan population yields a smaller satelUte-central offset. 

This is not exactly the same as the satellite population in the halo occu- 
pation models constrained as described above. However, in the HOD mod- 
els, the satellite mass as a function of infall mass is by definition preserved at 
smaller radii, and comparing the systems constrained by observed clustering 
at small radii < Ryii (l^t alone <C lOOkpc), which are near-merger, they are 
almost all in the orphan population in these particular SAMs. Regardless, the 
selection from the SAMs does not qualitatively change our conclusions. 
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Fig. 12. — Possible effects of allowing the satellite HOD to differ from the central HOD by a significant amount. Top Left: Average satellite versus average 
central galaxy mass (total baryonic mass, not stellar mass) at fixed maximum or infall halo mass (for the centrals, just Mhaio, for the satellites, their maximum 
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major merger rate (mergers per galaxy per Gyr, for galaxies with mass > M„ai), from convolving halo-halo merger rates with the above HOD for both centrals 
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low masses), owing to satellite-specific physics that lower 
Mgai (Mhalo) for satellites relative to its value at the time of in- 
fall/accretion. Semi-analytic models tend to predict that low- 
mass satellites (especially those about to merge - the popu- 
lation of interest) have substantially lower baryonic masses 
than central galaxies with the same halo mass at the time of 
the satellite accretion ("infall mass"). The baryonic mass dif- 
ference is about a factor of ^ 2 — 3 at the lowest masses, where 
disks are gas-dominated. We consider the difference in the 
predicted distribution (Mgai(Mhaio)) for satellites and centrals 

Note that this distinction does not necessarily appear in the stellar mass 
as well (S. White, private communication), but we focus on baryonic masses 
here since it more app ropri ately focuses on the relevant dynamical major 
merger events (but see § l5.1l below). 



at higher redshifts {z — I and z = 2) as well, and find that the 
gap between centrals and satellites, in the models, extends to 
higher masses at higher redshifts. 

Given this, we can estimate how large an effect the differ- 
ent satellite-versus-central distributions will have on merger 
rates. As above, we can compare the galaxy-galaxy merger 
rate on an even footing by adopting a fixed halo-halo merger 
rate (here that in Fakhouri & Ma 2008), and populating each 
halo according to the Mgai(Mhaio) relation. Unlike in the pre- 
vious section, however, we do not assume that both galaxies 
obey the same Mgai(Mhaio) relation. Rather, the primary and 
secondary subhalos separately follow the relations for central 
and satellite galaxies, shown here. 

Recall that differences between semi-analytic and semi- 
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empirical estimates of the Mgai(Mhaio) distribution for central 
galaxies have relatively little effect on merger rates. Thus, 
the variations in merger rates reflected in Figure [12] almost 
entirely owe to the difference (or lack thereof) between the 
Mgai(Mhaio) of satellites and that of centrals. This is in con- 
trast to the HOD-based model in which the Mgai(Mhaio) is 
nearly the same for centrals and satellites, leading to the same 
curves predicted in § 14.11 In the SAM, the median Mgai (Mhaio) 
can be significantly lower for satellites at the same infall/pre- 
accretion mass (especially at low masses). Thus even for a 
merger of exactly equal-mass halos (which, on average, must 
have equal-mass galaxies just before the halo-halo merger, i.e. 
at infall), the median expectation in the SAM is that the "sec- 
ondary" (whichever happens to have slightly less initial mass) 
will be only some fraction of the mass of the primary by the 
time the actual galaxy-galaxy merger occurs, if it were ini- 
tially gas-rich. The number of major mergers can therefore be 
substantially suppressed. 

It is clear how a larger discrepancy between Mgai(Mhaio) of 
satellites versus centrals at low masses leads to a stronger sup- 
pression of the merger rate. Likewise, at high redshifts, if 
the gap between central and satellite masses grows, the ma- 
jor merger rate will be further suppressed. At high masses, 
though, and in regimes where galaxies are gas-poor, both 
SAM and HOD constraints find a similar result where the 
A^s;ai(A^haio) distributions of satellites and centrals are nearly 
identical - as a result, their predicted merger rates are not dra- 
matically different. 

The results from the iBower et al.l (l2006l) and ICroton et al.l 
semi-analytic models are qualitatively similar. In con- 
trast, the^enii^analytic models of Somer villeetal. (2008^) 
and iKhochfar & Silkl (|20()9|) predict little discrepancy be- 
tween the Mgai(Mhaio) relations of satellites and centrals. This 
is because they do not include any mechanism for the satel- 
lite galaxies to lose a significant mass in gas before their final 
merger; we discuss this further below. 

4.2.2. Effective Mass Loss in Satellites: What Drives These 
Differences ? 

Satellite-specific physics in the semi-analytic models are 
clearly important for the predicted merger rates at low masses. 
Why do some models predict different Mgai(Mhaio) distribu- 
tions for satellites versus central galaxies? 

In general, the predicted differences in Mgai(Mhaio) between 
satellites and centrals relate to a well-known tendency for 
semi-analytic models to predict a satellite population that is 
"over-quenched." Satellites in these models tend to be too red, 
too gas-poor, too low-mass, and insufficiently star-forming 
relative to observed satellite galaxies. Figure [13] briefly il- 
lustrates this effect: we show first the observed fraction of 
passively evolving ("quenched") central and satellite galax- 
ies, as a function of galaxy stellar mass and halo mass. We 
compare these observations to the predicted passive fractions 
in semi-analytic models. 

The observational results are taken from 'Kimm et al.' 
(j2009), from their analysis of the Yang et al. (2005b, 2008) 
SDSS (z = 0) galaxy catalogs, with the passive population 
identified by the observed bimodal distribution of specific star 
formation rates at each mass. The observations show that, at 
fixed stellar and halo mass, the fraction of satellite galaxies 
that have "quenched" is not much different from that of cen- 
tral galaxies. At low stellar masses, in particular, most satel- 
lites are still star-forming (even those satellites in high-mass 
halos). Similar results have been found by other analyses 



and in a wide range of differe nt surveys, over a range of red- 
shifts from z = 0—l (see e.g.lGerke et alJl2007l:lH aines et al] 
12007'; 'Weinmann et al.' '2006a"b', '2009a'; 'Wang et al. 200'! 
Yang et al. 2008; Guo et al. 2009; van den Bosch et al. 2008). 
This agrees well with the expectation from semi-empirical 
models and HOD constraints (independently determined via 
the observed small-scale clustering of galaxies), that (at least 
massive) satellites should have properties that reflect their 
halo when it was last a central, but otherwise differ little from 
central galaxies. 

This is also reflected in galaxy gas content, shown in Fig- 
ure [T3] We show the observed gas fraction distribution of 
central and satellite galaxies (separately), as a function of 
stellar mass. We specifically show compilations of observed 
disk gas fractions as a function o f their stellar mass from 
iBell & de Jong ( l2000l) . iMcGaughl (l2005h . and iKannappanI 
(^2004^ . These observatio ns span both satel lite and central 
galaxies. The sample of IKannappanI (|2004|) . being simply 
mass-selected, is dominated by central galaxies at all masses. 
But from the other samples we can separate both central 
(chiefly massive field) galaxies and satell ite galaxies , in par - 
ticular the Ursa Major Cluster sample of iTuUv et al.l (1 19961) . 
This is ideal for our comparison, as it represents a rich, but 
not extremely massive (and thus unusual) group with a group 
velocity dispersion of ^ ISOkms^' -this is typical of the host 
halos of central/primary galaxies withMgai ^ 10'" — 10" M©. 
We stress that we are not trying to compare with galaxies in 
e.g. massive clusters, where the gas-rich galaxies tend to be 
much less massive than their cluster host - stripping and other 
processes may be efficient in such galaxies, but those would 
not represent major mergers in any model. This is also seen 
in many of the color or SFR - based studi e s above; for exam- 
ple, 'Wein mann et"an (l2009al ) : [Gug e t al.1 (i2009h: IWang et"d] 
(2009); van den Bosch et alJ (l2008h : iKang & van den BoschI 
(|2008j) find that only in the most extreme (rare) environments 
could there be significant gas depletion in satellites - over a 
wide range of satellite and halo masses, they find that there 
is at most very minor depletion (for major-merger candidates) 
and, in particular, that the amount of depletion does not in- 
crease or decrease with halo mass in any significant manner. 
This is true even in the central regions of groups, i.e. systems 
"about to merge," for which the predicted satellite-central gas 
mass difference in some models is maximized (see references 
above). The important result observed is, for galaxies with 
similar in similar-mass halos, which have some possibil- 
ity of being major merger pairs, there is not much difference 
in gas fractions between central and satellite systems. 

We compare both of these diagnostics - the "quenched" 
fraction and gas mass fractions in central and satellite galaxies 
- to the predictions of the semi-analytic models. The differ- 
ences are immediately apparent. The SAMs predict that es- 
sentially all satellites are quenched, especially those in halos 
above the critical "quenching mass" of ~ IO'^Mq. The av- 
erage predicted gas fraction of satellites at low mass is much 
less than that of equivalent central galaxies, in contrast to the 
observations. 

For clarity, we show in Figure [T3] just the pre- 
dictions fro r n the Ide Lucia & Blaizol (2007) model, but 
iKimm et"an (l2009t) and others have demonstrated that this 
satellite "overquenching" tendency is quite genera l in re - 
cent SAMs, includ i ng e.g. those of ICroton et alj ( 2006 ), 
'Monaco et al. (2007). Bower et al. (2006). Ka ng et alJ(l2005l) . 
Cattaneo et a l. (2006), and Bertone et al. (2007d! 

As noted above, we are re-stating a well-known result; the 
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Fig. 13. — Illustrations of the effect seen in Figure [T2I where low-mass satellites in many SAMs have suppressed baryonic masses. Left: Fraction of satellite 
and central galaxies that are star-forming (at z. = 0), as a bivariate function of galaxy stellar and halo mass. Top Left: Observed distributions from Kimm et aQ 
©09). Low -mass satellites have large star-forming fractions, similar to central galaxies. Bottom Left: Prediction from the de Lucia & Blaizot (2007) model. 
Satellites tend to be uniformly quenched. Right: Average gas fractions (and scatter in gas fractions) versus galaxy mass, for satellites and centrals. Top Right: 
Observed distributions at z = for field centrals I Kannappan 2004), and satellites in "normal" groups (not massive clusters, where satellite properties are 
significantly different; Bell & de Jong 2000). Barring extreme (rare) environments or very minor mass ratios, satellites are not strongly gas-depleted. Bottom 
Right: Prediction from the SAM for satellites and centrals. Model satelhtes are strongly gas-depleted. These are general trends in a number of semi-analytic 
models. 



same physical deficiencies are reflected in a number of other 
measures: in semi-analytic models, satellite mass functions 
tend to be under-massive, the fraction of blue satellites is 
lower than that observed, and the observed small-scale clus- 
tering of star-forming or blue galax ies tends to be poorly re- 
produced ( see references above and Cooper et al. 2006, 2007; 
Park et al l 1200 7: Blanton & Berlind 2007; Li et al. 2006; 
Coil et al.ll2008h . 

What, in these models, drives the "overquenching" of satel- 
lites? In most, once a system becomes a satellite, its halo gas 
reservoir is immediately added to that of its "parent" halo, and 
(as a consequence) no new cooling occurs onto the secondary. 
But if this removal of un-cooled, pristine halo gas were the 
only effect, the consequences would not be so pronounced. In 
all of these models, strong stellar feedback is assumed to af- 
fect star-forming galaxies - generally some mass loss of order 
several times the star formation rate is ejected from the disk, 
per unit stars formed. Given such strong mass ejection, in 
order to maintain observed (high) star formation rates in disk 
galaxies, the ejected gas must be recycled rapidly - either pro- 
cessed through the halo or some (for numerical reasons) des- 
ignated "reservoir" and returned quickly to the galaxy. But 
when a system is a "satellite" in these models (even if the 
mass ratio is close to 1:1), the mass is simply lost to the par- 



ent halo. The gas mass is therefore quickly depleted and the 
satellite quenches as it has ejected most of its gas.^"^ 

This can be a very large effect. A system with ~ 50% gas 
will lose ~ 1 /2 its baryonic mass in just a short time after 
accretion. Thus an initial 1 :2 major merger will become a 1 :4 
minor merger (or a more minor merger, if the parent "steals" 
this ejected mass). As a consequence, there will be almost no 
very gas- rich major mergers. 

This problem is further exacerbated in many of these mod- 
els by the assumption of instantaneous recycling from stellar 
evolution. When some stars form in the new satellite, the mass 
loss of ~ 50% that will occur over their lifetime is immedi- 
ately turned into gas, which can then suffer the consequences 
above. A more gradual mass loss over actual stellar lifetimes 
would, obviously, yield somewhat less rapid an evolution (es- 
pecially for major mergers, where the merger timescale might 
be only a Gyr). 

Since the consequences of this gas loss will clearly be 
stronger for systems with higher gas content, low-mass galax- 
ies are preferentially affected. This gives rise to the mass- 
dependent trend in the discrepancy between Mgai(Mhaio) pre- 

This is despite the fact that many of the models do not nominally include 
an additional explicit "ram-pressure stripping prescription" for the cold disk 
gas. 
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dieted for satellite and central galaxies by such models (shown 
in Figure fT2T l. where low-mass satellites have their masses 
more suppressed relative to centrals. Likewise, since at high 
redshifts all galaxies are more gas-rich, the discrepancy grows 
in the models with redshift. This, in turn, suppresses predicted 
merger rates at low masses and high redshifts. 

Differences between different semi-analytic models largely 
depend on exactly what prescription is used for stellar feed- 
back (controlling how efficiently the system self-quenches) 
as well as e.g. when and where this occurs and how it is 
connected to the lifetimes/stripping/identification of subha- 
los/satellites. For example, in some mo dels based on cos- 
mological dark ma tter merger trees (e.g. Ide Lucia & Blaizol 
l2007t iBower et al.1 12006). the system is labeled a "satel- 
lite" when the halos are first linked by the friends-of-friends 
(FOF) halo finder algorithm, which occurs for major pairs 
at anywhere from ^ 2 — 5/?vir- In such a case, these ef- 
fects can "strip" the secondary even before it enters the pri- 
mary halo. In contrast, simulations and observations find 
that, if any stripping occurs, it often only does so in the cen- 
tral, most dense and high-pre ssure regions of massive ha- 
los (Weinmannetal. 2006a b : IWang et al.ll2007t lYang et al.l 
2009c Weinmann et al. 2009?. 'van den B osch et al.1 120081: 
Kanp& van den Bosch 2008; McCarthy et alj|2008|). In fact, 
cosmological hydrodynamic simulations show that satellite 
galaxies tend to keep growing and accreting new material 
(both onto their subhalo and into the galaxy) until they are 
in the final sta ges of orbiting a relaxed massive system near 
its center (Si mha et al.l 120091) . Observations likewise see no 
significant difference between satellites in the outer regions 
of even massive groups, relative to centrals in similar-mass 
subhalos. The tendency of satellites to continue accreting 
in simulations is stronger at higher redshifts, and observa- 
tional constraints from z = — 2 require that the gas-rich, star- 
forming satellite fraction remain large over this entire inter- 
val (Tinker & Wetzel 2010), enhancing the discrepancy with 
these models at high redshift. 

Unsurprisingly, at high masses, where gas fractions tend to 
be small ('^ 10 — 20%), these effects are minimized. Thus the 
Mgai(Mhaio) distributions for satellites and centrals are similar, 
and there is little difference in predicted merger rates at these 
masses. 

4.2.3. Towards Resolution: What Physics for Satellites? 

It is important to recall that the vast majority of merg- 
ers by both number and mass density are in the field (es- 
sentially single halo-halo mergers) or small/loose group en- 
vironments (analogous to e.g. the Local Group). The ob- 
servations indicate that satellite-specific physics appear to 
be important only for systems with the smallest mass ra- 
tios /i„ai < 0.01 — 0.05 in these environments (specifically see 
e.g. IHaines et alJf2007l: IWeinmann et alJl200 9a: Kim m etaP 
lIooC and then only in the final stages and/or close pas- 
sages, not just at the momen t of crossing i^yir (refer- 
ences above and 'Yang et al."2009 Hvan"den Bosch et al.ll2008l: 
|Kang& van den Bosch 2008). Even if stripping were ex- 
tremely efficient in dense environments (e.g. for clusters, or 
for dwarf galaxies in massive halos such as the Local Group), 
such that it suppressed all such major mergers, it would have 
almost no effect on the overall merger rate. 

In contrast, much of the intuition used to design the pre- 
scriptions above is built up from satellite-specific studies 
based on massive clusters with dense, high-pressure ICM 
gas. But such systems contain only a small fraction of 



the bulge mass density of the Un iverse (< 5%; se e e.g. 
Carollo et al. 1998; Korm endv & Ke nnicutt 2004; Alle n et alJ 
i2006>;,Ball et al.,.2006,; .Driver et alJl2007uGadottii,2009i) . 

Indeed, with the possible exception of dwarf galaxi es in 
massive halos (D'Onghia et al.' 2009; Wang et alJl2009 '). the 
observations favor a "gentle" treatment of satellites. Consid- 
ering a volume-limited sample of satellites (i.e. one not dom- 
inated by extreme environments) that might become major 
mergers (i.e. are not so different in mass from their centrals), 
the observed specific star formation rates, colors, and gas frac- 
tions can be well-fitted (to first order) by a model in which 
there is no removal of the hot or cold satellite gas (including 
the "ejecta" or recycled material from stellar winds and su- 
pernovae). The only constraint is that the satellite halo stops 
growing, so e.g. the total baryonic mass supply is "frozen in" 
at the value it had when the subhalo was accreted. 

Such a simple model leads to a color and passive/quenched 
galaxy distribution as a bivariate function of galaxy mass and 
halo mass similar to that of centrals (with just slightly more 
"quenching" in low mass satellites), in better agreeme nt with 
observations (Weinmann et al. 2006a,b; va n den Bosch et alj 
120081; iKang & van den Bo sch 2008; iKimm et alj|2009l) . and 
similar gas fractions (as th e subhalo re tains eje cted stellar 
wind f eedback), as obser ved dBell & deT ong 2000; lMcGaughl 
120051; iKannappanI |2004|) . This also naturally leads to the 
observed result that satellite quenched fractions are a much 
stronger function of the satellite stellar mass, and weaker 
function of e.g. the local density or halo mass, than cur - 
rendy predicted (IWeinmann et al.ll2006bt iHaines et"ani2007h 
at bo th z = and high redshifts z ^ 1 — 1.5 ("Gerke et alj 
l2007h . And the observed weak dep endence of re d satellite 
fraction on distance from halo center ( IWeinmann et al.l2006ah 
is similarly obtained. The observationally favored gradual 
decline of satellite star formation rates relative to those that 
such systems had when they were central galaxies (see e.g. 
Wang et al. 2007) arises naturally in an effectively "closed- 
box" model, giver i the observed st ar formation rate-gas den- 
sity scaling laws dKennicuttl 1 1 998h . as systems deplete their 
original baryon reserve without new halo growth. 

High-resolution hydrodynamic simulations appear to give 
consistent results. Gas-rich systems are in the "cold accre- 
tion" regime: their halos are dense, with a cooling time shorter 
than the local dynamical time, so gas outflows may be effi- 
ciently halted by the dense ISM and more gas in new accretion 
is brought into the galaxy in cold filaments on the dynamical 
time. Especially in this regime, these cold flows dominate ac- 
cretion - the idea of spherical accretion from a quasi-static 
halo accounts for only a small amount of residual cooling in 
the most massive halos (essentially, this is the old "cooling 
flow" problem in massive clusters; but it contributes negligi- 
bly to the global cold gas budget or star formation rate den- 
sity). These cold filaments are sufficiently dense to survive 
disruption by the diffuse parent halo gas when systems be- 
come satellites; the simulations find that even relatively small 
satellites continue to hold onto their halo gas. 

In fact, the simulated satellites continue to accrete new 
material, and grow by cooling, well into their parent halo 
virial radius and for ^Gyr after become satellites (see e.g. 
ISimha et alj|2009l) . In this regime, satellite growth in the sim- 
ulations is indistinguishable from central galaxy growth, for 
the same halo/infall mass. At high redshift, the simulations 
suggest that even the "frozen in" models may quench satel- 
lites too severely; subhalos in the simulations continue to ac- 
crete new halo gas at a rate equal to the central galaxy well 
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into the halo. In effect, it is as if the system does not "know" 
it is a satellite until at very small radii from the primary center , 
just before the final galaxy-galaxy merger (ISimhaetal.1120091: 
Ikeres et al. 2Q09b a). 

Recall also that we are ultimately most interested in the 
case of major mergers: in these cases, it is clear that the dis- 
tinction of "satellite" versus "central" is somewhat semantic. 
Physically, the two are mutually gravitating (one is not mov- 
ing in the halo or potential of the other; rather the two mu- 
tually coalesce in a single halo/potential) and (by definition) 
have similar halo properties (and therefore similar halo gas 
masses, virial temperatures, etc.). Simple physical considera- 
tions (and simulations) imply that there should be no signifi- 
cant difference in the cooling/mass loss of one system relative 
to the other (let alone ram-pressure stripping of one and not 
the other). 

Of course, it is easier to point out these differences than to 
design physical prescriptions that effectively model them - let 
alone prescriptions that are robust across the entire regime of 
interest (from dwarf galaxies and/or small satellites in massive 
clusters to equal-mass pairs in field or Local Group-analogue 
environments). And as noted above, it is not enough to simply 
retain pristine halo gas associated with the satellite - ejected 
gas from stellar feedback should be recycled in some manner 
appropriate to how it is implemented in the different models. 

Consid erable progress is being mad e in this area . For 
example, iFont et al.r(l2008l) updates the iBower et all (l2006h 
model, with the primary change being a more detailed satel- 
lite cooling model. Specifically, the models are otherwise 
identical, but instead of removing subhalo gas instantly (as in 
[Bower et al. 2006), the revised model strips the gas from the 
assumed extended sub halo gas halo according to the model in 
[McCarthy etaH (I2008,) . as a function of e.g. the local gas den- 
sity (itself a function of radius from the primary, for assumed 
isothermal sphere gas profiles). There are still some sources 
of uncertainty in e.g. the exact wind treatment and efficiency 
of subhalo gas stripping (and the prescription adopted may not 
be applicable if most accretion comes via cold flows), but the 
agreement with the star-forming fractions of satellite galaxies 
is substantially improved. 

Figure [14] compares the merger rates predicted by t his re- 
vised model to those from the previous iBower et al.l (l2006h 
iteration of the model (with strong satellite quenching). We 
also compare the rates determined from HOD-based ap- 
proaches, and observational constraints. At high masses 
where systems are gas-poor and not star-forming, all the 
models (and observations) agree well. At low masses, 
where galaxies are gas-rich, however, there is an order-of- 
magni tude difference in the predicted merger rates. The re- 
vised IFont et al.l (l2008h yields merger rates more similar to 
the HOD-based models, as expected, since the primary adjust- 
ment to the model is intended to bring it into better agreement 
with the observed HOD statistics used as input in the HOD- 
based models. This yields significantly improved agreement 
with observational estimates of the merger rates. There are 
clearly remaining differences; however, they are at the factor 

2 — 3 level, typical as we have shown of the other differ- 
ences between the models. 

This is reinforced by comparison of Figures |2] & [12] Be- 
sides the iFont et al.l ( 1200 8) model discussed above, the two 
semi-analytic models that predict merger rates in best agree- 
ment with the halo occupation mo dels, over the e ntire galaxy 
mass range, are those from Khochfar & Silkl ( 2009 ) and 
ISomerville et al.l (l2008ah . These models do not include any 



mechanism for cold gas (including that affected by stellar 
fee dback) to be renioved f rom the satellite before a merger; 
the lSomerville et al] (l2008a ) model allows for some stripping 
of hot gas, but only inside of the parent virial radius where 
the timescale for a major merger is short and so there is little 
effect on the satellite mass. In Figure[2] we showed that these 
models are precisely those which predict little or no discrep- 
ancy between the Mgai(Mhaio) relation of satellites and cen- 
trals. This does not mean that there are not other, perhaps 
very significant, differences between the models (nor does it 
necessarily mean that the same answer is being obtained for 
the correct reasons); but it does imply that, controlling for dif- 
ferences in predicted satellite properties removes a significant 
factor driving different predictions for merger rates. 

5. THE OTHER SIDE: DEFINITIONS OF MASS AND EFFECTS OF 
MERGERS ON MORPHOLOGY 

5.1. The Definition of Mass Ratio 

Throughout this paper, we have focused on the "merger 
rate," with respect to mass ratios defined in a specific man- 
ner Observationally, of course, some of these are more or 
less accessible (the galaxy-galaxy luminosity or stellar mass 
ratios being most so, the total dark matter halo-halo mass ra- 
tios least so). But physically, the same "number" of mergers 
will have very different implications for galaxy growth, mor- 
phologies, bulge formation, and star formation, depending on 
the mass ratio involved - and the merger rate is not the same 
for different definitions of the mass ratio. 

Figure [TS] illustrates some of these caveats. First, consider 
the merger rate itself, within the same model, but given differ- 
ent definitions of merger mass ratio. There are three obvious 
choices of mass ratio, commonly used: galaxy-galaxy stellar 
mass ratio (/i* = M*,2/A/*, i), galaxy-galaxy baryonic (stel- 
lar-Hcold gas) mass ratio (/igai = Mgai,2/A^gai, i), and halo-halo 
mass ratio (^haio — A^haio.a/A^haio, i)- We have shown how the 
shape of the Mgai — Mhaio relation means that the number of 
major mergers in terms of /^gai can be quite different from the 
number of mergers in terms of /ihaio- But since galaxy gas 
fractions are not uniformly constant, /igai and /i* will also not 
be the same. 

In Figure [15] we plot the galaxy major merger rate (major 
mergers per galaxy per Gyr) at z = 0, as a function of galaxy 
stellar m ass. We do so for a single model - the "default" 
model in iHopkins et al.l (l20To|) using the halo-halo me rger 
rates from'Fakhouri & Ma! (|2008|);|F^ouri et al.' (2010") and 
abundance-matching HOD method from Conrov & Wechsle^ 
( |2009|) . But the results are qualitatively the same regardless 
of model. We show the major merger rate for different defini- 
tions of a "major" merger: a merger with halo-halo mass ratio 
Mhaio > 1/3. with galaxy stellar-stellar mass ratio /i* > 1/3, 
and with galaxy baryonic-baryonic mass ratio //gai > 1/3. 
The halo-halo rate is only weakly dependent on mas s, a re - 
sult well-known from cosmological simulations (see § 13.2.11 ). 
The stellar-stellar rate, on the other hand, depends strongly on 
mass - this owes to the shape of M*(Mhaio)- At low masses, 
galaxy stellar mass is a steep function of halo mass, suppress- 
ing major mergers {M^, cx M^aio' ^ ^'^ halo-halo merger is 
just a 1:9 stellar-stellar merger); at high masses, stellar mass 
is a shallow function of halo mass, enhancing major merg- 
(M, oc mV 



ers 



'halo' ^° ^ halo-halo merger is a 1:3 stellar- 



stellar merger). The effect is very large - an order of mag- 
nitude at the extremes of the mass function - and cannot be 
neglected. The baryonic-baryonic rate is quite similar to the 
stellar-stellar rate at high masses (not surprising, since high- 
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Fig. 14. — Comparison of the predicted major galaxy merger rates (as Figure[2) in two similar semi-analytic models, 'Bower et al.'(200fl) and'Font et alJKOOS'). 
The SAMs are otherwise identical, but Fon t et al. (2008) implement a more physical, less severe satellite stripping/quenching model, improving the agreement 
with e.g. the observations in Figures [l2ll3l We compare the predictions to estimates directly adopting the observational HOD constraints, and to observed merger 
rates. Rates at h igh M^-^i, where cooling/gas fractions are negligible, are unchanged. Rates at low masses are significantly higher in the revised model, as expected 
from Figure nUand agree more closely with direct HOD models. 
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Fig . 15. — Left: Importance of the definition of "mass ratio" for merger rates. We show the major (p > 1/3) merger rate per galaxy at z = 0, for galaxies 
at a given stellar mass, in terms of the merger halo-halo mass ratio, galaxy stellar-stellar mass ratio, baryonic-baryonic (stellar-Hcold disk gas) mass ratio, and 
"dynamical" mass ratio (gravitational mass, i.e. baryons plus the tightly bound, non-stripped dark matter inside the baryonic radii; see Equation|4). The different 
ratios behave very differently as a function of mass. Halo mass ratios are most easily predicted, but very different from galaxy mass ratios. The baryonic or 
dynamical mass ratio is most relevant for e.g. morphological disturbances and bulge formation, but difficult to directly observe. Stellar mass ratios are accessible 
in e.g. pair samples, but under-estimate the rate of "damaging" major mergers at low masses by factors ~ 3. Right: Impact of mergers in high-resolution 
hydrodynamic simulations from Hopkins et al. (2009a). We show the remnant bulge-to-total mass ratio after a merger with mass ratio /x; to remove ambiguity 
in the merger mass ratio and rate, all systems have just one merger, with the systems initialized so that fi = fj,f = /ihalo = Mgal = Mdyn- We show simulations 
with different orbital parameters (symbol types), gas fractions (color, from /gas = in red to /gas = 1 in black), and structural parameters (variation in disk sizes, 
initial B/T, and halo concentrations in solid). Other parameter variations include surveys of stellar wind models (magenta), black hole feedback (violet), and gas 
equations of state (grey). At fixed, perfectly-defined /i, there is factor of 2 — 3 variation in the amount of bulge formed by mergers owing to these parameters - 
neglecting their full distributions will inherently limit any prediction of the "impact" of mergers. Merger impact is also cleai'ly continuous across the major/minor 
merger (p = 1/3) distinction. 



mass galaxies are uniformly gas-poor). At low masses, how- 
ever, it lies somewhere between the stellar-stellar and halo- 
halo rate; this is because low-mass galaxies have increas- 
ing cold disk gas masses, even while their stellar mass falls 
rapidly with halo mass. Thus the total baryonic mass is not 
as steep a function of halo mass as is the stellar mass, and 
major mergers are less strongly suppressed (at low masses, 
Mgai oc M^^^^^'^ , so e.g. a 1:2 halo -halo merger remains a 
major baryonic-baryonic merger; see lAvila-Reese et ani2008l : 



iBell et alJ l2003al: iMcGaughl l200l IStaAetan l2Q09h . The 
differences are again very non-trivial - at low masses, the 
baryonic-baryonic merger rate is a factor of 3 larger than the 
stellar-stellar merger rate. Fo r furth er discussion of these dis- 
tinctions, we refer to lStewarl (l2009l) . 

These various ratios not only relate differently to observable 
quantities; they also have a number of physical relationships 
to the impact of galaxy mergers. Using halo-halo mergers 
as a proxy for galaxy-galaxy mergers is clearly not a good 
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approximation in either models or observations, except just 
at ^ L* where the mass function turns over Physically, the 
halo-halo mass ratio is also not the most important quantity 
for e.g. bulge formation. This is because, by the time the final 
merger/galaxy coalescence takes place, most of the secondary 
halo mass will already have been stripped and/or mixed with 
the primary halo mass. In general, the outer parts of the halo 
(which contain much of the mass), being extended and low- 
density relative to the baryonic galaxy, are not able to induce 
significant perturbations or torques on that galaxy. 

The stellar-stellar mass ratio has the advantage of being 
observable, as well as identifying at least a portion of the 
tightly bound material in the galaxies that can actually in- 
duce structural perturbations and bulge formation. However, 
at low masses or high redshifts, where gas fractions are non- 
negligible, stellar mass alone clearly misses a very important 
quantity - the amount of cold, rotationally supported gas in 
galaxy disks. Observations indicate that in these regimes the 
cold gas is gravitationally dominant over the stars in a large 
fraction of the galaxy population. From a gravitational per- 
spective, in terms of the perturbation induced, it makes no dif- 
ference whether the same merging disk is 100% gas or 100% 
stars (although the end effects will be d ifferent and gas may b e 
distributed differently from stars; see iHopkins et alj|2009al) . 
even though one would be a stellar-stellar major merger and 
the other would not appear in any observed merger catalog. 

Thus, attempts to estimate the "impact" of mergers (e.g. the 
amount of bulge they might form), or even the total number 
of mergers a galaxy experiences, based on observed merger 
fractions identified with stellar mass or luminosity cuts, will 
under-estimate the true impact of mergers by a factor of at 
least ^ 3 at low masses (M* < IO'^Mq). More directly sam- 
pling the total baryonic content of these galaxies is necessary 
to overcome this limitation; in the meantime, empirical re- 
sults should be corrected for this distinction or at least rec- 
ognize its very large importance at low masses. And gas 
fraction is critical for understanding the consequences of a 
given merger; the dynamics of mergers are fundamentally dif- 
ferent and bulge formation less efficient in gas-dominated, 
as opposed to stellar-dominated, systems (see IHopkins et al] 
|^09a). At high masses, > 10" M©, at least at low redshifts 
z < 1, such a procedure does not produce a result very differ- 
ent from the baryonic-baryonic merger rate. 

But even the baryonic mass ratio is not really representa- 
tive of the gravitational mass that induces perturbations and 
forms bulges, heats disks, and scatters stars in mergers. What 
is actually desired is the surviving, tightly bound total gravi- 
tational mass of material that is involved in the last couple of 
passages in the final merger. In high-resolution simulations of 
a large nu mber of mergers, spanning a large parameter space, 
ICox et ai] (f2008) find that this can be roughly approximated 
by the sum of the galaxy stellar mass, the cold disk mass, and 
the dark matter mass that is tightly bound to the baryons and 
thus will be able to resist stripping - approximately that in- 
side a radius 2 — 4Re (where Re is the effective radius of 
the baryonic galaxy) or ~ 1 rj (where = Rm/c is the NFW 
scale radius of the halo, with c the halo concentration). Note 
that since this dark matter retention at small radii is largely a 
baryonic effect, it will not appear in collisionless simulations. 

We therefore define a "dynamical" mass Mdyn, and corre- 
sponding mass ratio fidyn- 

Mdyn +M„^,+Mdm{< 3 Re) 

^dyn =Mdyn,2/Mdyn,l . (4) 



This is a significantly better approximation to what is most 
important for bulge formation and the dynamical properties 
of galaxies than any of the mass ratios discussed above. Fig- 
ure[T5]shows the corresponding major merger rate, as a func- 
tion of primary stellar mass, where a major merger is de- 
fined by fi^yn > 1 /3 (to calculate the dark matter content in- 
side 3 Re, we assume the halos follow NFW profiles with a 
concentration-halo mass relation from iBullock et al.l (1200 ll) . 
and that the gala xy Re follows the relation as a function of 
stellar mass from 'Sh en et al.l 12003)). Roughly, the result is 
similar to that for baryonic mass ratio, but it is slightly closer 
to the halo-halo merger rate at both high and low masses (un- 
surprisingly, since it includes baryons plus a halo contribu- 
tion). 

Of course, this mass is difficult to extract observationally. 
However, it is readily accessible in most semi-empirical and 
semi-analytic models. Adopting such a definition for mass 
ratio, as opposed to the simpler definition of baryonic mass 
ratio, is important for predictions at low and high masses. In 
fact, many such models still adopt either just the stellar-stellar 
mass ratio or the halo-halo mass ratio as their proxy for the 
"impact" of a merger - this is not only clearly not the most 
physical quantity, but as Figure [15] illustrates, it can lead to 
an order-of-magnitude over or under-estimate of the conse- 
quences of mergers. 

5.2. The Importance of Other Parameters 

Even with perfect knowledge of the merger rate and dis- 
tribution of mass ratios, however, we stress that the im- 
pact of mergers cannot be determined at better than fac- 
tor ~ 2 accuracy, absent a large amount of additional in- 
formation. Figure [15] illustrates this. We consider a large 
suite of high-resolution, hydrodynamic galaxy merger sim- 
ulations. They include realistic progenitor galaxy models 
(gasH-stellar disk-nbulge-nhalo), cooling, star formation, and 
feedb ack from star formation and BH growth in a multi-phas e 
ISM (ISpringel & HernquistI l2003al: iDi Matteo etalllioOSh 
These a re presented and discussed in detail in IHopkins et al.l 
( I2009ai) ; for now, we use them to illustrate a couple of sim- 
ple points. Although it may not be the most representative for 
galaxies at extreme masses, we wish to remove the already- 
discussed ambiguities in merger rates and mass ratio defini- 
tions. Therefore, every galaxy undergoes exactly one iso- 
lated merger with the simulations all initialized such that the 
"mass ratio" is the same by any of the definitions above, i.e. 
A* = A'haio = A** = Mgai — Mdyn- We run the simulations un- 
til they are fully relaxed, and then quantify the bulge-to-total 
mass ratio B/T in the remnant stars - the property most of- 
ten desired in relation to mergers. We consider mass ratios 
fi = Q.l — I, and vary a large number of other parameters in 
the simulations, including galaxy gas fractions, orbital param- 
eters, feedback prescriptions, and structural properties of the 
initial disks, halos, and bulges. 

First, it is very clear that there is no special "division" at 
the "major" merger threshold of /i = 1/3. In many models, 
it is simply assumed that major mergers completely destroy 
disks, while minor mergers leave them intact. Likewise, in 
many observational studies, simple assumptions such as this 
are used to relate the number of major mergers to some "ex- 
pected" number of bulges. But in fact, to lowest order, bulge 
formation is continuous, with the fraction of the initial disk 
destroyed (remnant B/T) scaling sls B/T ^ u- This is dis- 
cussed in greater detail in Hopkins et al. | (|2009al). with similar 
results found in a number of studies (iBournaud et al.l 120051: 
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[Johansson et al.ll2009bl: lYounger et al.ll2008h. and sonie of th e 

consequences are discussed in lHopkins 

Such a simple assumption regarding the efficiency of mergers 

can lead to systematic factors of ^ 2 — 10 differences in the 

total bulge mass formed (depending on e.g. the galaxy mass), 

and clearly skews the distribution of mergers impacting bulge 

formation. 

Also clear in Figure [15] is a large scatter in remnant B/T 
at fixed despite the fact that the mass ratio and number of 
mergers are perfectly well-known and defined for these sim- 
ulations. This stems from the other parameters varied in the 
simulations. Most important are the orbital parameters and 
gas fractions. We consider a wide variety of different orbital 
parameters (denoted by different symbols). Because stellar 
scattering and the gaseous angular momentum loss that drives 
dissipation, starbursts, and the formation of the bulge from 
the galaxy gas are dominated by r esonant processes, orbital 
param eters have a huge impact ( Barnes & Hernquis3ll99lL 
Il996l) . It is well-known that prograde mergers tend to be 
much more destructive (form much larger bulges) than ret- 
rograde mergers, and we see that in Figure [15] We also see a 
continuum of beh avior in between, as has been noted in many 
numerical studies ('H ernauistll989t |Hernauist & Mihosl ll995 



Naab & Burkert 20031: ICox etllMO OS; Younger et al 

i\ 



2008 



Bournaud et al.il2005riHopkins et al.li2009ai Johansson et alj 
2009blah . The effect is large, with e.g. a typical 1:3 merger 



leading to B/T k, 0.5 in the most prograde mergers and 
just B/T ^ . 25 in c orresponding retrograde mergers. In 
iHopkins et all (l2009ah . it is noted that gas loses angular mo- 
mentum and falls to the galaxy center, contributing to a nu- 
clear starburst and bulge mass, within an initial critical radius 
R (xl/{l —ip cos 6*), where tp ^ 0.5 — 0.6 depends on struc- 
tural properties and 9 is the initial disk inclination (0 — is 
prograde, = tt retrograde) - i.e. angular momentum loss is 
efficient out to a radius w 3 — 4 times larger in prograde cases, 
relative to retrograde cases. 

Likewise, the gas fraction has a similar effect - the most 
gas-poor mergers lead to very efficient bulge formation while 
the most gas-rich mergers lead to very little bulge formation 
even in extremely major mergers. Figure [15] shows a large 
number of /i = 0.5 — 0.8 cases with only B/T ^ 0.2 — 0.3, 
provided that the gas fractions at the time of the merger are 
high ( /gas ~ 0.6 — 0.8). Again, this has been seen in a number 



of works iSpringel & Hernauisil2005 


;lRobertson et al.ll2006t 


Governato et al. 2009: Richard et al. 


l2010c Johansson et al.l 


2009blat). In Hockins et al. (2009al), these simulations are 



used to show that, to lowest order, the efficiency of gas an- 
gular momentum loss scales with a factor oc (1 — /gas) - i-e. 
in a merger with 50% gas, the ultimate gas angular momen- 
tum loss per unit mass is (on average) half as efficient as in a 
merger with just 10% gas. 

Other parameters lead to further scatter in B/T at fixed 
/i. Varying galaxy structural p arameters, for exampl e, can 
stabilize or destabilize systems (ISpringel et al.ll2005al) . The 
presence of initial bulges in a certain mass and scale radius 
range suppr esses gas inflows and disk heating on first pas- 
sages (Mihos & Hernquist 1994, 1996). Allowing gaseous 
disks to be more extended than stellar disks leads to less effi- 
cient angular momentum loss in the gas (because that angular 
momentum loss is usually dominated by transfer from the gas 
to the stars). Changing the prescriptions for star formation, 
stellar feedback, pressurization of the multi-phase interstellar 
medium, and cooling will lead to more or less efficient gas 
compression and star formation in the earlier phases of the 



merger - this changes the absolute amount and spatial dis- 
tribution of the gas at the time of the final merger, and thus 
the efficiency with which gas is channeled into a central star- 
burst and bulge. In sufficiently minor mergers, tidal stripping 
of the external portions of the halo can effect some of the 
gas and stars, and this can be accelerated by ram-pressure ef- 
fects if we model a halo with a very dense, high-pressure hot 
gas background (analogous to e.g. massive clusters). There 
is even the possibility, especially in very low-/i mergers, of 
complete tidal disruption of the secondary. This may be im- 
portant for the overall mass budget in the most massive galax- 
ies (although that depends on the precise separation between 
'galaxy' and intra-cl uster or intra-g roup light, which is poorly 
defined; see e.g. Dol ag et al.ll2009l) . 

It is not our purpose here to address these effects in de- 
tail; we refer to Hopkins et al. (2009a) for a more thorough 
discussion. However, we wish to stress that if any physical 
consequence of mergers is desired, then even perfect knowl- 
edge of the merger rate and mass ratio distribution is only 
sufficient for an order-of-magnitude estimate of the impact of 
mergers. Further knowledge of e.g. the distribution of merger 
orbital parameters and gas fractions is necessary to improve 
estimates of the impact of mergers at the factor ~ 2 — 3 level. 
And to improve to an accuracy much beyond a factor ^ 2, 
knowledge of the much more detailed structural properties of 
galaxies and highly uncertain effects of feedback would be 
required. 

6. SUMMARY AND DISCUSSION: THE UNCERTAINTY BUDGET 

A simple comparison of different predictions of the galaxy- 
galaxy merger rate from e.g. different semi-analytic models 
and simulations demonstrate that the predictions vary by an 
order of magnitude (Figure [2]). We have attempted to survey 
the sources of uncertainty and systematic differences between 
various theoretical attempts to predict the galaxy merger rate, 
to identify what drives these differences and address how 
progress can be made. 

6.1. Dark Matter and Dynamics 

All models depend similarly on the "background" dark mat- 
ter merger rate. However, with proper caution in adopting 
definitions, and modern convergence in high-resolution cos- 
mological simulations, the uncertainties in this rate can be 
reduced to the factor ~ 2 — 3 level. In other words, if how 
galaxies populate halos is appropriately fixed, there is rela- 
tively little uncertainty in the merger rate owing to the dy- 
namics of the dark matter (Figure [TJ. The relevant quantities 
considered here include: 

Cosmology: Nominally, changing the cosmology within 
the stat istical uncertainty in modern observational constraints 
(iKoma tsu et al. 2009) leads to factor 1.5 differences in the 
merger rate. However, most of this owes to changes in the 
mass function which can be normalized out - if e.g. the galaxy 
population is normalized so as to match the observed stellar 
mass function and large-scale bias, then the resulting differ- 
ences in the merger rate are negligible compared to the other 
sources of uncertainty below (Figure[T]i. 

Halo-Halo Merger Rate Determinations: Again, nomi- 
nally different halo-halo merger rate determinations differ by 
factors of ^ 2. Some of this stems from e.g. the inherent am- 
biguity in masses and extents of halos. However, much owes 
to the definitions adopted when attempting to fit/quantify the 
instantaneous merger rate "function." If these definitions are 
properly accounted for, or if halo merger trees are used di- 
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rectly to track galaxies (rather than defining halo merger rates 
at all), the real uncertainties are small, a factor ^1.5. Effects 
of baryons contribute similarly small 10 — 20%) uncertain- 
ties in the overall halo-halo merger rate (Figure[3]i. However, 
care is needed with definitions - defining halo mergers simply 
in terms of "instantaneous" mass ratios, and/or applying cuts 
which can be useful for constructing galaxy merger trees, or 
using timestep-sensitive mass ratio definitions can lead to an 
apparent (artificial) suppression of major halo-halo mergers 
by an order-of-magnitude or more (Figures|4l|5]). This requires 
careful consideration in simulation-based semi-analytic mod- 
els. 

Subhalo versus Halo Merger Rates: Defining the merger 
rate not when halos first merge into e.g. a larger friends-of- 
friends group, but when subhalos merge/are destroyed into 
the central (primary) group subhalo, may be more represen- 
tative of galaxy-galaxy mergers. Doing so, however, intro- 
duces additional uncertainties owing to e.g. how subhalos are 
identified. It is also important, in this case, that the subhalo 
"mass" (for merger mass ratio purposes) be defined as the "in- 
fall" or maximum pre-accretion/pre-stripping mass (i.e. max- 
imum mass before the system became a subhalo), otherwise 
the mass of all systems at merger, by definition. Differ- 
ent calculations from various simulations, with different sub- 
halo identification/destruction criteria, yield rates converged 
to within a factor ^ 2 (Figure|6ll. Using hydrodynamic simu- 
lations to "tag" mergers (using the galaxies therein as the ideal 
subhalo/halo "tracers") yields consistent results. 

Merger Timescales: In models without subhalo resolution 
or with Press-Schechter merger trees, galaxy-galaxy merg- 
ers are often assumed to be "delayed" with respect to the 
halo-halo merger by a timescale given by e.g. the dynami- 
cal friction time (this approximates the subhalo evolution). 
We show that calibrations of such timescales, from high- 
resolution galaxy merger simulations, are relatively well- 
converged (Figure |7ll. Adopting one of these calibrations, 
we show that the implied merger rate agrees well with that 
obtained from full tracking of subhalos, with similar fac- 
tor ^ 2 uncertainties (Figure [8]l. For major mergers, both 
methods agree well with the simple halo-halo merger rate - 
this is because the timescale for such a merger to complete 
(^ fHubbie) is short Compared to the time between such merg- 
ers (-- fHubble)- 

However, we show that adopting an artificially long 
merger timescale, or adopting e.g. the simple Chandrasekhar 
timescale with a normalization higher than the specific numer- 
ical calibrations here, can suppress high-redshift mergers by 
factors ^5 — 10. We also caution that these calibrations are 
designed for the time from halo-halo merger - they are not 
constructed for application in "hybrid" models, where subha- 
los are followed in simulations down to some resolution limit, 
and then a merger time is applied to the residual "orphan" 
based on its instantaneous radius and mass. We show that in 
such cases, if the system is lost to resolution at large radius 
> Q.lRyii, the calibrated formulae are not self-consistent, and 
can lead to over-estimates of the total merger timescale by 
factors of - 2 - 8 (Figure|9ll. 

6.2. Baryonic Physics 

Controlling for some of the caveats and definitions above, 
the dominant uncertainties in predicted merger rates owe to 
variations in how galaxies populate halos. The differences can 
be identified with the shape of the assumed halo occupation 
distribution (HOD) - essentially, the distribution of galaxy 



masses in given (sub)halo masses, Mgai(Mhaio)- Some model 
for this is necessary to translate halo-halo mergers (or halo- 
subhalo mergers) into galaxy-galaxy mergers. We therefore 
consider how these physics are modeled and how they lead to 
variations in the merger rate, in three classes of models with 
very different approaches towards modeling the galaxy-halo 
connection: 

Semi-Empirical (HOD-based) Models: In semi-empirical 
models the HOD is adopted explicitly from observational con- 
straints. As such, it is subject to the attendant uncertainties 
and limited to the dynamic range observed. At low redshifts 
and galaxy masses within factors of several around ^ L,, the 
constraints are tight and different methods agree well - result- 
ing uncertainties in the galaxy-galaxy merger rate (holding 
halo-halo merger rates fixed) are a factor 1 .5 (Figures [T] & 
[Tol l. The uncertainties grow to a factor ~ 2 at the lowest and 
highest masses, and z ^ 1 — 2. Above z ~ 2, the uncertainties 
grow very rapidly: there are not sufficient observational con- 
straints on the HOD to make strong statements about galaxy 
merger rates. Predicted merger rates and pair counts from 
such models agree well with direct observations over a stellar 
mass range ^ 10^' ^ — lO^'^Mf?^ and redshifts z = — 2 (se e 
Figures above andlS tewart et al.ll2009al: iHopkins et al.ll20I0h . 
At this level, it appears, there is no tension between ACDM 
merger rates and merger/pair counts. To the extent that other 
models disagree significantly with the observed merger frac- 
tions, it should owe to issues in the model baryonic physics 
leading to a Mgai(Mhaio) distribution different from that ob- 
served. 

Cosmological Hydrodynamic Simulations: In hydrody- 
namic simulations, the Mgai(Mhaio) distribution is predicted 
in an a priori manner based on the cooling, star formation, 
and feedback models implemented in the simulation. Unfor- 
tunately, it is not yet possible to run large-volume cosmologi- 
cal hydrodynamic simulations (needed to quantify the galaxy- 
galaxy merger rate) with the spatial and mass resolution and 
detailed prescriptions for feedback from stars and black holes 
that it is becoming clear are necessary to form "realistic" 
galaxies. Simulations available that do not include feedback 
yield poor agreement with the observed HOD and stellar mass 
function, predicting a relationship closer to the efficient star 
formation limit (Mgai = //,Mhaio)- Relative to what the merger 
rates would be, taking the same dynamics and merger lo- 
cations but re-populating galaxies with masses chosen to fit 
the observed HOD (stellar mass and clustering), these simu- 
lations tend to over-predict merger rates at low masses and 
high redshifts, and under-predict rates at high masses and 
low redshifts by factors ^3 — 5, as well as over-predicting 
the relative importance of minor versus major mergers at all 
masses (FigurefTOb. Because it is the shape of the Mgai — Mhaio 
relation that is most important, simply re-normalizing pre- 
dicted masses by a uniform factor will not correct for these 
effects. Re-normalizing all masses to their "correct" masses 
given some observed HOD is an improvement, but care is still 
needed, since the incorrect masses and morphologies will af- 
fect quantities such as the dynamical friction time. 

Semi-Analytic Models: In semi-analytic models, the diffi- 
culties and expense of simulations are replaced by use of ana- 
lytic prescriptions, given some background dark matter popu- 
lation, to predict ultimate galaxy properties. Such models are 
adjusted to give good agreement with the galaxy stellar mass 
function (and clustering) at z = 0; as such, some agreement 
with the Mgai(Mhaio) distribution of central galaxies (which 
dominate the stellar mass function at all masses) is implic- 
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itly guaranteed. Indeed, we find that adopting the predicted 
SAM HODs for central galaxies, instead of the empirically 
determined HOD, yields no systematic difference and scatter 
within the same factor ~ 2 allowed by different observational 
constraints (Figures [TOll. At high redshifts, the uncertainties 
in both grow. Some SAMs yield growing discrepancies rela- 
tive to semi-empirical models, directly related to issues such 
as e.g. the known tendency of SAMs to over-predict the abun- 
dance of low-mass galaxies at high redshift, but these are still 
within the factor 3 — 5 level at z < 3 (Figure [TTTi. 

However, the HOD of satellite galaxies in semi-analytic 
models is extremely sensitive to prescriptions for cool- 
ing, stellar feedback, and halo mixing/stripping in satellites. 
Moreover, satellite masses are not strongly constrained by the 
stellar mass function, so there is no implicit guarantee/check 
that these are correct, and more detailed comparison with e.g. 
observed group catalogs and small-scale clustering must be 
used to calibrate the models. 

In detail, it is well-known that most SAMs have difficulty 
reproducing the observed properties of satellite galaxy popu- 
lations: they tend to predict satellite galaxies that are under- 
massive and over-quen ched, relative to observa tio ns (see Fig- 
ures [H] & [T3] and e.g. I\yeinmann etani2006allbt IWang et all 
|2007; Kim metal.|[2009r and references therein). This can 
occur even in models where the small-scale clustering of 
satellites is ove r-predicted, relative to observations (see e.g. 
iGuo et ai]l2010b . who find both effects, although the cluster- 
ing discrepancy is probably due to the large value of ag used in 
the simulation). In most of the models, satellites lose their en- 
tire halo gas reservoir at the moment they become such - even 
in major mergers (and even when the moment of "becoming" 
a satellite occurs at several times the primary virial radius). 
Moreover, the combination of simple stellar wind feedback 
and cooling models often leads to the satellites also losing 
almost all of their cold/disk gas reservoir 

As such, at low masses where gas fractions are important, 
initially major and even equal-mass mergers can easily be- 
come minor mergers in the model, potentially suppressing the 
predicted merger rate by a large factor Models which yield 
more "over-quenched" satellite populations tend to yield cor- 
respondingly smaller merger rates (Figures [12] & \14[ . Cor- 
recting for these differences, for example by enforcing that 
satellite galaxies obey a similar HOD to central galaxies, 
or by adjusting the prescriptions for cooling onto satellite 
galaxies such that they better reproduce the observed color 
and star formation rate distributions of satellites, leads to 
larger merger rates that converge with the merger rates pre- 
dicted from semi-empirical models (Figure O. Whether or 
not this is necessary for matching e.g. close pair counts and 
small-scale clustering is unclear (see e.g. the comparison in 
liCitzbichler & White 2008; Guo et al. 2010b), but it demon- 
strates an important uncertainty in predictions of merger rates. 

6.3. Impact of Mergers on Galaxies 

Of course, even with perfect knowledge of the merger rate 
and distribution of mass ratios under some definition, it is not 
trivial to relate this to either observable or physical quantities 
such as merger fractions or the amount of bulge mass formed 
in mergers. We consider how two basic aspects of this relate 
to uncertainties in merger rates and their consequences: 

Mass Ratio Definitions: As outlined in Stewart (2009), a 
halo-halo major merger is not necessarily a galaxy-galaxy ma- 
jor merger, and vice versa. But there are also several means of 
defining galaxy-galaxy mass ratio, including e.g. the stellar- 



stellar mass ratio /i, and baryon (stellar mass plus cold gas 
within the galaxy disk)-baryon mass ratio /Zgai. At high (low) 
masses, the merger rate in terms of stellar-stellar mass ratio 
is enhanced (suppressed) by an order of magnitude relative 
to the halo-halo merger rate (Figure [TSll. In comparing e.g. 
model predictions to stellar mass ratio-selected pair samples, 
clearly the stellar mass ratio is most applicable. However, 
in comparing to luminosity-ratio selected pair samples, or to 
morphological studies, where the total baryonic mass is what 
matters, the merger rate in terms of the baryon-baryon mass 
ratio is more relevant. At high masses, this behaves similarly 
to the merger rate in terms of the stellar-stellar mass ratio. At 
low masses, however, galaxies are increasingly gas rich - thus 
many mergers with minor stellar-stellar mass ratios in fact 
have major baryon-baryon mass ratios, and the baryon-baryon 
major merger rate is a factor ^ 3 higher than the stellar-stellar 
merger rate. This may, in fact, explain at least part of the 
well-established fact that morphology-inferred merger rates at 
low masses tend to be system atically higher than pair coun t- 
inferred merger rates (see e.g. iLopez-Sanjuan et al]l2009allR 
and references therein). 

Although often a superior definition, even the baryonic 
mass ratio misses the tightly-bound dark matter that is im- 
portant for the dynamics of final merger, which will not be 
stripped because it is within the baryonic radii. We there- 
fore propose a new "dynamical" mass and mass ratio which 
approximates the most important quantity in high-resolution 
simulations, and includes both baryons and some dark mat- 
ter Merger rates in terms of this quantity behave similarly 
to baryonic major mergers, but with somewhat (^ 20 — 50%) 
higher (lower) rates at low (high) masses (Figure [TSll. 

If one wishes to infer the amount of bulge or other "dam- 
age" done by major mergers but uses either the stellar-stellar 
mass ratio or halo-halo mass ratio, the estimate can easily be 
systematically incorrect at the factor ^3 — 10 level. This is 
critical for observational studies seeking to infer the amount 
of bulge formed by majo r mergers, using a s tellar-mass-ratio 
selected sample (see e.g. iBundv et al.ll2009l) : at low masses, 
there might be ^ 3 times as many "damaging" mergers as 
given by this statistic. Also, many analytic models simply use 
the stellar (or even halo) mass ratio as a criterion for determin- 
ing the impact of a merger - the systematic errors introduced 
by this assumption can be larger than those from halo mis- 
identification, complete ignorance of subhalos, assignment of 
incorrect merger timescales, or discrepancies between model 
and observed stellar mass functions. Yet despite the consid- 
erable literature on those problems, there has been relatively 
little focus on the adoption of better mass ratio proxies. 

Otlier Parameters and Merger 'Impact': Even when 
controlling for the above uncertainties, we show that at per- 
fectly well-defined merger rates and merger mass ratios, there 
is a large variation in the physical effects of mergers, owing 
to other parameters such as the merger dynamics and orbit, 
merger gas fractions, and initial structural properties of the 
merging systems (Figure [Tsll. At the same merger mass ratio, 
a prograde orbit can build twice as much bulge as a retrograde 
orbit, and will lead to much more dramatic tidal features and 
distortions observable in morphology-selected samples. Also, 
the merger timescale will be significantly different at moder- 
ate/small radii < lOOkpc, so pair-selected samples will also 
see biased distributions. A gas-poor major merger will typi- 
cally violently relax the entire stellar disk and funnel the gas 
entirely into a starburst that builds central bulge mass, but a 
very gas-rich merger will experience very inefficient angular 



30 



Hopkins et al. 



momentum transfer, suppressing the burst mass by a factor 
^ (1 — /gas), a factor ~ 3 — 5 at low masses and high red- 
shifts. More subtle properties lead to scatter and offsets at a 
smaller, but non-negligible level. Many of these are discussed 
in detail in Hop kins et alJ ('2009a). 

From a purely empirical perspective, without knowledge of 
these properties, it is difficult to assess the impact of merg- 
ers (in terms of e.g. the bulge mass formed) at better than the 
factor '^2 — 3 level. From the theoretical perspective, neglect- 
ing quantities such as orbital parameters and, especially, gas 
fractions, in forward-modeling merger remnants again will in- 
troduce systematic uncertainties that are larger than any of the 
uncertainties from modeling the dark matter distribution, sub- 
halos, merg er timescales, and the like. Many of these are out- 
lined in Hopkin s et al.l(l20 09b); at high masses, the systematic 
uncertainties are less severe, but at low masses, where systems 
tend to be gas-rich, they can be factors ^ 5 — 10 in the total 
bulge mass formed. 

Another important point is that mergers are continuous - 
there is no special division at the traditional 1:3 distinction 
between "major" and "minor" mergers. In many models, 
and in many observational assessments of merger "effects," 
it is assumed that major mergers are completely destructive, 
while minor mergers do no damage. This simple assumption 
introduces systematic factor ^ 2 errors in the total merger 
"damage budget" an d total bulge formation efficiency (see 
iHopkins et alJl2()Tol) : the effect in terms of skewing which 
mergers do "more" or "less" for bulge formation can obvi- 
ously be severe. 

6.4. Observational Constraints and Outlook 

Testing these models and determining the true merger rate 
in a robust manner will of course ultimately depend on ob- 
servations. Continued observations of the merger rate, bear- 
ing in mind the caveats above, are of obvious importance. 
In improving such estimates, calibration of specific samples 



to high-resolution A^-body simulations, specifically with mock 
observations matched to the exa ct selection and m ethodology 
adopted, will be critical (see e.g. lLotz et alj|2008ab . And even 
with such a calibration, the discussion above makes it clear 
why the relation between observed merger rates and bulge 
buildup (depending on a number of secondary properties not 
directly observed), in detail, must ht forward-modeled. 

Tighter observational constraints on the halo occupation 
distribution - from e.g. group catalogs, kinematics, weak lens- 
ing, and clustering - in particular at low masses and at high 
redshifts, will directly improve the semi-empirical models, 
and put strong constraints on the a priori galaxy formation 
models in the areas that have greatest effect on predicted 
merger rates. As we have shown, constraints on satellite pop- 
ulations specifically will be valuable. However, the satellite 
populations of particular interest are not the historically well- 
studied extreme cases of e.g. dwarfs in the local group or 
low-mass galaxies in Virgo and massive clusters (from which 
much of our intuition regarding dynamical friction, stripping, 
and satellite gas exhaustion comes). The case of interest for 
merger rates is that of major mergers (i.e. near equal-mass 
galaxies) in field or loose group environments (i.e. systems 
analogous to the local group, but where either the Milky Way 
or Andromeda is the "satellite" of interest). Also, at low stel- 
lar masses and high redshifts, large uncertainties remain in 
galaxy gas masses, and these matter as much or more relative 
to the stellar mass in assessing the "impact" of mergers. 
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